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A METHOD FOR DETERMINING THE MOD-/ BEHAVIOUR OF 

RECURSIVE SEQUENCES 


C. KRATTENTHALERt AND T. W. MULLER* * 


Abstract. We present a method for obtaining congruences modulo powers of a prime 
number p for combinatorial sequences whose generating function satisfies an algebraic 
differential equation. This method generalises the one by Kauers and the authors 
[Electron. J. Combin. 18(2) (2012), Art. P37] from p = 2 to arbitrary primes. Our 
applications include congruences for numbers of non-crossing graphs and numbers of 
Kreweras walks modulo powers of 3, as well as congruences for FuB-Catalan numbers 
and blossom tree numbers modulo powers of arbitrary primes. 


1. Introduction 


The purpose of the present article is to generalise the basic theory concerning con¬ 
gruences for combinatorial sequences developed in [T2] for 2-power moduli to arbitrary 
prime power moduli, and to discuss a number of further applications of our method. 
Given a prime number p, consider the (formal) power series 



( 1 . 1 ) 


n>0 


The series •hp/) is easily seen to be transcendental over Z; this can be shown, for 
instance, by generalising (in a straightforward manner) the density argument given in 
m Lemma 1], In Lemma CD below, we provide a different, and more elegant, argument, 
which however presupposes a number of non-trivial facts concerning certain auxiliary 
series Hb l j )2 ,....hr{ z ) established in Section [3] 

While <& p (z) is transcendental over Z, it is easily seen to be algebraic when considered 
over finite rings Z/p 7 Z, and the first focus of our paper is on polynomial identities for 
$ p (z) modulo a given p-power, which are of minimal degree ; cf. the next section. A (still 
stubbornly unproven) conjecture concerns this minimal degree of $ p (z) modulo p 7 , see 
Conjecture @] in the next section. As in [121 Sec. 2] for the case of the prime p = 2, 
we are only able to establish that the conjectured minimal degree constitutes a lower 
bound, see Lemma CD 

The basic idea of our method for producing and establishing congruences modulo 
p-powers for combinatorial sequences is to express solutions of appropriate (formal) 
differential equations (in technical terms: differentially algebraic power series over the 
integers satisfying a uniqueness condition) as polynomials in & p (z), with coefficients 
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that are Laurent polynomials in z. For each such differential equation, our method may 
be implemented to (semi-)automatically yield such polynomial expressions in $> p (z) for 
arbitrary p-power moduli, provided such expressions do in fact exist. 

However, such a procedure as described in the last paragraph would be without any 
practical value, if we could not, at the same time, provide an efficient algorithm for 
extracting coefficients from powers of $ p (z). This is the topic of Section [3] Briefly, the 
solution to this problem we propose is to (i) expand powers of & p (z) into Z-linear sums 
of certain auxiliary series 

H bl , b2 ,..., br (z) := ^+^ 2 +-+^ ; ( 1 . 2 ) 

ni>n2>--->n,->0 


(ii) to show that the set consisting of 1 and the series (11. 2D with parameters bi, b 2 , ..., b r 
relatively prime to p are linearly independent over the ring Z [z, z^ 1 ], as well as over the 
rings (Z/p 7 Z)[z, z -1 ]; (iii) to prove that a series H blb2 _. :br (z), without restrictions on its 
parameters b±, b 2 , ..., b r , can always be expressed as a linear combination over Z [z, z _1 ] 
of these linear independent series; and (iv) to demonstrate that coefficient extraction 
from a series (11.21) with all parameters bi,b 2 ,..., b r relatively prime to p can be done 
effectively ; cf. Equation (13.ip . Corollary 0 Proposition [9J and Remark [8] 

A detailed description of our method is given in Section [4], while the remaining 
sections treat a number of applications beyond the case where p = 2: in Section [51 we 
obtain congruences modulo arbitrary 3-powers for the number of non-crossing graphs 
with a given number of vertices; Section 12 deals with congruences modulo 3-powers 
for the number of Kreweras walks of given length in the plane; Section [7] treats the 
Fufi-Catalan numbers 


F(n; k) : = 


n 


kn 

n — 1 


modulo powers of a prime number p, provided that the parameter k is itself a power 
of p: finally, the last section studies the number 


B{n\p) 


p + 1 ( pn \ 

n((p — 1 )n + 2) \n — l) 


of p-ary blossom trees with n white nodes modulo arbitrary powers of p. See, in par¬ 
ticular, Theorems EU [H HSJ M (H EDI E2 and ESI 

We would like to point out that Rowland and Yassawi HZ] have developed a com¬ 
pletely different method for proving congruences. Their method applies to diagonals of 
rational functions, and is based on the construction of automata. Our method is com¬ 
plementary: there is a considerable intersection of applications (in particular, if one 
also adds the variation [13] of the method presented here), but neither can our method 
cover all diagonals of rational functions nor can their method cover all differentially 
algebraic series. Applications of our method not covered by P2J can be found in [ 121 
Sections 8-14], 


2. The p-power series <3> p (z) 

The series $ p (z) in (II.ip is the principal character in the method for determining 
congruences of recursive sequences modulo p-powers which we describe in Section [4j For 
the sake of better readability, we shall usually suppress the index p and simply write <f>(z) 
for < F p (^). It is not hard to see that this series is transcendental over Z[z] (this follows 
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for instance from Lemma Q] below). However, if the coefficients of <E>(,s) are considered 
modulo a p-power p 7 , then < h(^) obeys a polynomial relation with coefficients that are 
polynomials in z. The focus of this section is on what may be said concerning such 
polynomial relations, and, in particular, about those of minimal length. In the proofs 
of Lemma |T] and Proposition [2] below, we shall already make use of some fundamental 
results which will only be established later in Section [3] This will not create any circular 
arguments since Lemma [Tj and Proposition [2] will not be needed in Section [3j 

Here and in the sequel, given power series (or Laurent series) f(z) and g(z ), we write 

f(z) = g(z) modulo p 7 

to mean that the coefficients of z l in f(z) and g(z) agree modulo p 7 for all i. 

We say that a polynomial A(z, t) in z and t is minimal for the modulus p 7 , if it is monic 
(as a polynomial in t), has integral coefficients, satisfies A(z, <E>(z)) = 0 modulo p 7 , and 
there is no monic polynomial B(z,t) with integral coefficients of t-degree less than that 
of A(z,t) with B(z,<&(z)) = 0 modulo p 7 . (Minimal polynomials are not unique; see 
Remark 0) Furthermore, we let v p (a) denote the p-adic valuation of the integer a , that 
is, the maximal exponent e such that p e divides a. 

The lemma below provides a lower bound for the degree of a polynomial that is 
minimal for the modulus p 7 . 

Lemma 1 . If A(z,t) is minimal for the modulus p 1 , then the degree d of A(z,t) in t 
satisfies v p (d\) > 7 . In particular, the series <E>(z) is transcendental over’L[z\. 

Proof. Given positive integers b\, b 2 ,..., b r , define the series Hbub 2 ,...,b r ( z ) by 

H bl , b2 ,... A (z) := z b lP ^ + b^ 2+ ... +brp r,r _ 

ni>n 2 >--->n r >0 

By (13.ip and Proposition [9J a power $ d (z) can be expressed as a linear combination of 
series H bl ,b 2 ,...,b s (z) and 1 with coefficients in Z[z\. Furthermore, all bf s in these series 
are relatively prime to p, and b\ + b 2 + • • • + b s < d. In particular, we have 

$ d (z) = d\ H\ 1 i(z) + other terms, ( 2 . 1 ) 

with d repetitions of 1 in the index of H. Since, by Corollary [TJ the series H blt b 2 ,..., bs {z) 
and 1 are linearly independent over (Z/p 7 Z)[z], and (with d repetitions of 1) 

does not appear in the expansion of < 1 >e (z) for e < d, a monic polynomial of degree d in 
< F(^) can only vanish modulo p 7 if d\ = 0 (mod p 7 ), or, equivalently, if v p (d\) > 7 . This 
establishes the assertion. □ 

Proposition 2. Minimal polynomials for the moduli p, p 2 ,p 3 ,... ,p p+1 are 

t p — t + z modulo p, 

(t p — t + z) 2 modulo p 2 , 


(t p — t + z) p modulo pP, 

(t p — t + z) p — p p_1 (t 2 — t + z) modulo pP +1 , 
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Proof. By (13.111 . we have 


b\\b 2 \ - ■■ b \ Hbl,b2, - ,br ^ 

r= 1 bi,...,b r >l 

6lH-h b r =p 

= ~ r ‘*'+pit E .ua 

n> 0 r=2 bi,...,6 r >l i- r ' 


biH- \-b r =p 


or, equivalently, 


*(z)+z= P -£ ■£ .*■<*>• 

r=2 6 1 ,...,6 r >l 1 2 r 

blH-h b r =p 


( 2 . 2 ) 


Clearly, this implies the claim modulo p, and, by potentiation, as well the claims mod¬ 
ulo p l for i — 2, 3,... ,p. 

Taking both sides of (12.211 to the p-th power, we obtain 


($p(z) - <f>(E) + zf = pP 


v 

E 

r= 2 bi,...,b r >l 

b\-\ - \-br — p 


E 


(P ~ !) ! 

bi\ b 2 \ ■ ■■b r \ 


H p 

n b 1 ,b 2 ,...,b, 


(z) + O (f * 1 ), 


where the expression O (p p+1 ) subsumes terms all of whose coefficients are divisible 
by p p+1 . If, in the sum over r on the right-hand side, we use the congruences N p = 
N (mod p), which is valid for any integer N, and 


K.k .= .*.(*)+ 0 (p), 

then, in combination with ( 12 . 2[1 read from right to left, we obtain the claim on the 
minimal polynomial modulo pP +l . □ 


Remark 3. Minimal polynomials are highly non-unique: for example, the polynomial 

(t p -t + z) 2 +p(t p -t + z) 

is obviously also a minimal polynomial for the modulus p 2 . 

Based on the observations in Proposition [2] and Lemma [T] we propose the following 
conjecture. 


Conjecture 4. The degree of a minimal polynomial for the modulus p 1 , 7 > 1 , is the 
least d such that v p (d\) > 7. 


Remark 5. (1) Given the p-ary expansion of d, say 

d = d 0 + di ■ p + d 2 ■ p 2 + ■ ■ ■ + d r ■ p r , 0 < di < p — 1, 

by the well-known formula of Legendre [T5, p. 10], we have 


»(rfO=E v =E 


d 


i=\ L 

r i 


P 


1=1 


E<i 

,i=0 


iP 


i—t 




i =1 1=1 


J2 J2 d ^~ 

1=1 i=l 

p l — 1 d — s(d ) 

f " l p — 1 p — 1 


(2.3) 
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where s(d) denotes the sum of digits of d in its p-ary expansion. Consequently, an equiv¬ 
alent way of phrasing Conjecture 0 is to say that the degree of a minimal polynomial 
for the modulus p 7 is the least d with d — s(d) > (p — 1 ) 7 . 

(2) We claim that, in order to establish Conjecture HJ it suffices to prove the conjecture 
for 7 = ^~Ep, 5 = 1,2,.... If we take into account Lemma [1] plus the above remark, this 
means that it is sufficient to prove that, for each 5 > 1 , there is a polynomial As(z,t) 
of degree p s such that 

Ag(z, <3>(z)) = 0 modulo (2.4) 

For, arguing by induction, let us suppose that we have already constructed Ai(z,t), 
A 2 (z,t),..., A m (z, t ) satisfying (12.41) . Let a be a positive integer which is divisible by p 
and 

a = • p + a 2 ■ p 2 H- ha m - p m , 0 < a* < p - 1 , 

be its p-ary expansion. Then we have 

m m 

Jj4^(z,$(z)) =0 modulo =i ,(«--(«))/(P-i). ( 2 . 5 ) 

5=1 <5=1 

On the other hand, the degree of the left-hand side of (12.51) as a polynomial in <£>( 2 ) is 

E m S 

8=1 a SP = a - 

Let us put these observations together. In view of (12.31) . Lemma Q] says that the 
degree of a minimal polynomial for the modulus p 7 cannot be smaller than the least 
integer, dA 1 ' 1 say, for which d (r,> — s(d^) > (p — 1 ) 7 . (We point out that d^ must be 
automatically divisible by p.) If we take into account that the quantity a — s(a), as a 
function in a, is weakly monotone increasing in a, then ( 12 . 5 j) tells us that, as long as 

^( 7 ) < (p — l)p + (p — l)p 2 + • • • + (p — l)p m = p m+1 — p, 

we have found a monic polynomial of degree d^\ B^(z, t ) say, for which B 1 (z : < h(^)) = 0 
modulo p 7 , namely the left-hand side of (12.51) with a replaced by d^ 7 \ to wit 

m ^ 

B ^(z,t) =Y[A d 5 s (. z,t ), 

5=1 


where d^ = d^ ■ p + d^ • p 2 + • • • + d\n ■ p m is the p-ary expansion of d^\ Hence, it 
must necessarily be a minimal polynomial for the modulus p 7 . 

Since 

p m+1 — p — s(p m+1 — p) = p m+1 — p — (p — 1 )m, 

we have thus found minimal polynomials for all moduli p 7 with 7 < p +_~ p — m. Now 
we should note that the quantity a — s(a) makes a jump from p m+1 — p — (p — 1 )m to 
pm+i _ w j ien we m ove from a = p m+l — p to a = p m+1 (the reader should recall that it 
suffices to consider a’s which are divisible by p). If we take A^Zjt), which has degree 
p . p m = then, by (12.41) . we also have a minimal polynomial for the modulus 

(p(p m -i)/(p-i)) p = p (p m+ 1 -p)/(p-1) 


and, in view of the preceding remark, as well for all moduli p 7 with 7 between 


p m + l_ p 

P-1 


m + 1 and P m+1 p . 

p-i 
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So, indeed, the first modulus for which we do not have a minimal polynomial is the 
modulus ph’ m+1 “P)/h’~ 1 ) +1 = p(p m+1 -i)/(p-i). This is the role which A m+ i(z,t) (see (12.4p 
with m + 1 in place of 8) would have to play. 

The arguments above show at the same time that, supposing that we have already 
constructed Ai(z , t), A 2 (z,t),..., A m (z , t), the polynomial A^(z, t ) is a very close “ap¬ 
proximation” to the polynomial A m+ i(z,t) that we are actually looking for next, which 
is only “off’ by a factor of p. In practice, one can recursively compute polynomials 
Ag(z,t ) satisfying (12.41) by following the procedure outlined in the next-to-last para¬ 
graph before Lemma 0 in the next section. It is these computations (part of which are 
reported in Proposition [2]) which have led us to believe in the truth of Conjecture [4] 


3. Coefficient extraction from powers of $(z) 


In the next section we are going to describe a method for expressing formal power 
series satisfying certain differential equations, after the coefficients of the series have 
been reduced modulo p k , as polynomials in the p-power series <&(z) (which has been 
discussed in the previous section; for the definition see (II. 11)1 . the coefficients being 
Laurent polynomials in z. Such a method would be without value if we could not, at 
the same time, provide a procedure for extracting coefficients from powers of <3>(z). The 
description of such a procedure is the topic of this section. 

Clearly, a brute force expansion of a power <h A (z), where K is a given positive integer, 
yields 


K 


«%-) = E 


it! 


r =1 &i,...,6 r >l 

b±-\ -K b r =K 


bi\b 2 \ ■ ■ ■ b r \ 


II, 




(2). 


(3.1) 


where 





z b iP n l+b 2 p n 2+-+b r p nr 


ni>n 2 >--->n r >0 


The expansion (13. ip is not (yet) suited for our purpose, since, when bi,b 2 ,... ,b r vary 
over all possible choices such that their sum is K, the series H blb2 _ br {z) are not linearly 
independent over the ring Z[z, z~ x ] of Laurent polynomials in 0 over the integer^, and, 
secondly, coefficient extraction from a series H bljb2t ___ tbr (z) can be a hairy task if some of 
the fej’s are divisible by p. 

However, we shall show (see Corollary [7]) that, if we restrict to fej’s which are rela¬ 
tively prime to p, then the corresponding series H bl b2i ^, r (z), together with the (trivial) 
series 1, are linearly independent over Z[z,z~ 1 ], and there is an efficient algorithm to 
express an arbitrary series H bltb2t ... tbs (z), without any restriction on the feds, as a lin¬ 
ear combination over Z[z, of 1 and the former series (see Proposition [9]) . Since 
coefficient extraction from a series H blib2i ___^ T {z) with all b^s relatively prime to p is 
straightforward (see Remark [8]), this solves the problem of coefficient extraction from 
powers of $(z). 

As a side result, the procedure which we described in the previous paragraph, and 
which will be substantiated below, in combination with (13.11) . provides all the means 
for determining minimal polynomials in the sense of Section [21 as explained in Item (2) 
of Remark [7] at the end of that section, it suffices to find a minimal polynomial for the 


■^The same is true for an arbitrary commutative unital ring in place of the ring Z of integers. 
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modulus p( pS 1 \ 5 — 1, 2,.... For doing this, we would take a minimal polynomial 
A$_i(z, t ) for the modulus p( p WAp- 1 )^ expand Aff^z, t ), and replace each coefficient 
c at g of a monomial z a t 13 in A^^z, t ) by c a ^ + pO^Pl/Cp -1 )^^ where x a ^ is a variable, 
thereby obtaining a modified polynomial, Bg_i(z, t ) say. Now we would substitute $(2;) 
for t, so that we obtain Bg_ 1 (z, $(-7)). Here, we express powers of d>(z) in terms of the 
series H ai>a2> _ !ar (z) with all cq’s relatively prime to p, and collect terms. By reading the 
coefficients of z 1 H aua2t _ t a r (z) in this expansion of Bs-i(z,$(z)) and equating them to 
zero modulo p^ v ~ 1 F(p^ 1 ) ) we produce a system of linear equations modulo p( p W/Cp- 1 ) 
in the unknowns x aj p. By the definition of Ag~i(z, t), after division by p^ p this 

system reduces to a system modulo p, that is, to a linear system of equations over the 
held with p elements. A priori, this system need not have a solution, but experience 
seems to indicate that it always does; see Conjecture [4j 

We start with an auxiliary result pertaining to the uniqueness of representations of 
integers as sums of powers of p with multiplicities, tailor-made for application to the 
series H aua2t _ j(lr (z). 

Lemma 6. Let d,r,s be positive integers with r > s, c an integer with \c\ < d, and let 
a±, 02, ■ ■ ■, a r respectively b\, b 2 ,..., b s be two sequences of integers, none of the cq ’s or 
bi’s divisible by p, with 1 < a, < d for 1 < i < r, and 1 < bi < d for 1 < i < s. If 

a lP 2rd + a 2 p 2(r - 1)d + • • • + a r p 2d = b lP ni + b 2 p n2 + • • • + b s p n ° + c (3.2) 

for integers ni,n 2 , ■ ■ ■, n s with n\ > n 2 > ■ ■ ■ > n s > 0, then r = s, c = 0, = bi, and 

Hi = 2 d(r + 1 — i) for i = 1,2,... ,r. 

Proof. We use induction on r. 

First, let r — 1. Then s = 1 as well, and (13.2j) becomes 

a±p 2d = bip ni + c. (3.3) 

If n 1 > 2d, then the above equation implies 

aip 2d = c modulo p 2d+1 . (3.4) 

Consequently, the integer c must be divisible by p 2d . By assumption, we have |c| < d. 
Since d < p 2d , the only possibility is that c vanishes. But then it follows from (13.41) that 
Gq is divisible by p, a contradiction. 

If d < rii < 2d, then it follows from (13.31) that c must be divisible by p ni . Again 
by assumption, we have |c| < d < p d < p ni , so that c = 0. But then (13.31) cannot be 
satisfied since bi is assumed not to be divisible by p. 

If 0 < n\ < d, then we estimate 

bip ni + c<d (p d + l) < (p d - 1 )(p d + 1) < p 2d , 
which is again a contradiction to (13.31) . 

The only remaining possibility is ri\ = 2d. If this is substituted in (13.3j) and the 
resulting equation is combined with |c| < d < p 2d , then the conclusion is that the 
equation can only be satisfied if c = 0 and aq = b±, in accordance with the assertion of 
the lemma. 

We now perform the induction step. We assume that the assertion of the lemma is 
established for all r < R, and we want to show that this implies its validity for r = R. 








C. KRATTENTHALER AND T. W. MULLER 


Let t be maximal such that n t > 2d. Then reduction of (13.2p modulo p 2d yields 

b t+ ip nt+1 + b t+2 p nt+2 +-b b s p na + c = 0 modulo p 2d . (3.5) 

Let us write b ■ p 2d for the left-hand side in (13. 5[) . Then, by dividing (13.2p (with R 
instead of r) by p 2d , we obtain 

a 1 p 2{R - 1)d + a 2 p 2( ' R ~ 2)d H- Va R -ip 2d = +b 2 p n2 ~ 2d H-f b t p nt ~ 2d +b-a R . (3.6) 

We have 


0 < b < p~ 2d d [p 2d ~ l + p 2d ~ 2 -|- b p 2d ~ s+t + l) 

< p~ 2d d (j^{p U ~ p 2d ~ s+t ) + l) < d. 

Consequently, we also have \b — a R \ < d. This means that we are in a position to apply 
the induction hypothesis to (13.6p . The conclusion is that t — R— 1, b — a R — 0, a* = bi, 
and rii = 2 d(R + 1 — i) for i — 1, 2,..., R — 1. If this is used in (13. 2p with r = R, then 
we obtain 

anP 2d = c 
or 

a R p 2d = b R p nR + c, 

depending on whether s = R — 1 or s = R. The first case is absurd, since c < d < 
P 2d < dRP 2d - On the other hand, the second case has already been considered in (13.31) . 
and we have seen there that it follows that c = 0, a R = b R , and n R = 2d. 

This completes the proof of the lemma. □ 

The announced independence of the series Lfai,a 2 ,...,o r (^) with all ad s relatively prime 
to p is now an easy consequence. 

Corollary 7. For any commutative unital ring R , the series Lfai,a 2 ,...,a r (^) 5 with all ai’s 
relatively prime to p, together with the series 1 are linearly independent over R[z, z -1 ]. 
In particular, they are linearly independent over (Z/pZ)[z, z -1 ], over (Z/p 7 Z)[^, z _1 ] 
for an arbitrary positive integer 7 , and over 7L[z , z~ x \. 


Proof. Let us suppose that 


N 


q 0 {z) + J2^ H ad\a^,...,a^( Z ) = 0 ’ 


i= 1 


(3.7) 


where the qi(z) , s are non-zero Laurent polynomials in z over R , the rf s are positive 
integers, and a^, j — 1, 2,..., r*, i — 1, 2,..., N, are integers relatively prime to p. We 

may also assume that the tuples (oq', a 2 \ ..., ), i — 1, 2,..., N, are pairwise distinct. 

Choose io such that r* 0 is maximal among the rf s. Without loss of generality, we may 
assume that the coefficient of z° in q io (z) is non-zero (otherwise we could multiply both 
sides of (13.71) by an appropriate power of ^). Let d be the maximum of all a d> 's and 
the absolute values of exponents of z appearing in monomials with non-zero coefficient 
in the Laurent polynomials qi(z), i = 0,1 ,N. Then, according to Lemma El with 
r = r io , aj = Oj\ j = 1, 2,..., r io , the coefficient of 


y af 0 ' 1 p 2rd +af 0 \ 2( ' r ~ V>d -\ - \-ai°' > p 2d 
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is 1 in H (i 0 ) (i 0 ) (i 0 )(^), while it is zero in series z e H (; 0 ) (,: 0 ) (<o) (^), where e is a 
non-zero integer with \e\ < d, and in all other series z e H (») (*) «(<£), i = 1,..., 

“l ,u, 2 r i 

io — 1, io + 1, • • •, N, where e is a (not necessarily non-zero) integer with \e\ < d. This 
contradiction to (j3.7|) establishes the assertion of the corollary. □ 


Remark 8. Coefficient extraction from a series H aita2t ,__ tCLr (z) with all afs relatively prime 
to p is straightforward: if we want to know whether z M appears in H ai a2 ^ ^ ar (z), that 
is, whether we can represent M as 


M = aip ni + a 2 p n2 +-b a r p nr 

for some ni,ri 2 , ■ ■ ■ ,n r with rii > n 2 > • • • > n r > 0, then necessarily n r = v p (M), 
n r -1 = v p (M — a r p nr ), etc. The term z M appears in H ai!a2t _ !(lr (z) if, and only if, the 
above process terminates after exactly r steps. This means, that, with n r ,n r _i,... ,n\ 
constructed as above, we have 

M - ( a s p ns H-b a r -\p rir ~ 1 + a r p nr ) > 0 

for s > 1, and 

M — ( aip ni + • • • + a r _ip nr_1 + a r p nr ) = 0. 

It should be noted that, given a±, a 2 ,..., a r , this procedure of coefficient extraction 
needs at most 0(log M) operations, that is, its computational complexity is linear. 


It remains to show that an arbitrary series Ht n ,b 2 ,...,b s ( z ) can be expressed as a linear 
combination over Z [z^z^ 1 ] of the series 1 and the series H aitCl2 _ ^ ar (z), where all a^s are 
relatively prime to p. 


Proposition 9. For any positive integers bi,b 2 ,.. .,b r , the series i?b 1 ,6 2 ,...,6 s (0) can be 
expressed as a linear combination over Z[z\ of the series 1 and series of the form 
H ai ,a 2 ,...,a s {z)i where all ai’s are relatively prime to p, s < r, and a\ + a 2 + • • • + a s < 

b\ + b 2 + • • • + b r . 


Proof. Given H blb2t bs {z) where not all bf s are relatively prime to p, we let h be 
maximal such that b b = 0 (mod p). If there is no such h, then there is nothing to do. 
If there is, then we apply the appropriate case of 


H bl ,...,b h - 1 ,pb' h ,b h+ i,...,b r ( z ) 

Hbi,...,bh-i,b' h ,bh+i,...,b r (^) ~b H bl ,...,b h _i+b' h ,b h +i,..., br (*0 ■^6i,...,6h_i,p6| 1 +6h+i,...,6r (^)) 

if 1 < h < r, 

= < H b ^ b2: _ A (z) - H pb > i+b2} ' A (z), if 1 = h<r, 

H bi,., bh -i,v h (z) + H bi,..,b h -i+b' h ( z ) - z b ' h H bu ..., bh _ x {z), if 1 < h = r, 

k (z) — z 6 i, if 1 — h — r. 

(3.8) 


These identities are due to Hou mi, and they are straightforward to verify. The reduc¬ 
tion (13.8p is now applied to the series which arose during the first application of (13.81) . 
This process is repeated until all series H aitCL2 ^ aa {z) in the obtained expression have 
the property that each index a* is relatively prime to p. From (13.81) . it is obvious that 
all these series satisfy a x + a 2 + • • • + a s < b\ + b 2 + • • • + b r . This proves the assertions 
of the proposition. □ 
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To conclude this section, let us provide an illustration of the above discussion. We set 
ourselves the task of determining the coefficient of z 29 ‘ 398301914493 i n for $( 2 ) = 

£ n>0 , 3 ". In order to accomplish this task, we first express <f> 5 (^) in terms of series 
H ai ,...,a r ( z ) with all ads being relatively prime to 3. Indeed, by means of the expansion 
(13. lj) and the algorithm described in the proof of Proposition [9l we have 


® 5 (^) = 120Hi,i,i,i,i(z) + 60^2,1,1,1^) + 60^1,2,1,1^) + 60^1,1,2,1^) + 60ffi,i,i,2(^) 

+ 30^2,2,1 (^) + 30^2,1,2 ( 2 ) + 30 ^ 1 , 2 , 2 ( 2 ;) — 15i?4 i i(z) — 15Hi^[z) — 9H 5 (z) 

+ 6 OP /1 n(z) + 30H 2 1 (z) + 30Hi 9 (z) — 20zHi 1 (z) — 10zH 2 (z) + 10Hi(z) — lOz. 

(3.9) 

Now we have to answer the question, in which of the series H ait ... tar (z) that appear 
in this expansion of 4> 5 (,s) do we find the monomial ~ 297398301914493 . Using the algorithm 
described in Remark [SI we see that 

297398301914493 = 3 30 + 3 29 + 3 28 , 

= 2 • 3 29 + 2 • 3 29 + 3 28 , 

= 2 • 3 28 + 3 28 + 2 • 3 28 + 183014339639688, 

= 3 30 + 2 • 3 28 + 2 • 3 28 , 

= 4 • 3 29 + 3 28 , 

= 3 30 + 4 • 3 28 , 

= 3 29 + 3 29 + 2 • 3 29 + 3 28 , 

= 3 30 + 3 28 + 3 28 + 2 • 328, 

= 5• 3 28 + 183014339639688, 

= 2 ■ 3 29 + 3 28 + 137260754729766, 

= 3 28 + 2 ■ 28 + 228767924549610, 

= 2■ 3 28 + 251644717004571, 

= 3 28 + 274521509459532. 


Here, the first line shows that ^ 297398301914493 appears in i7i,i,i(z) (thereby making it 
impossible to appear in i(z), # 2 , 1 , 1 , 1 ( 2 ), or H 121 i{z)), the fifth line shows that 

it appears in H^i(z), and the sixth line shows that it appears in i7 4 ) 1 (z), while the 
remaining lines show that it does not appear in any other term in the expansion of 
4> 5 (z) displayed in (13.911 . Hence, by taking into account the coefficients with which 

the series -Hi, 1 , 1 ( 2 ), H^i(z), and Hi^{z) appear in this expansion, the coefficient of 
^,297398301914493 j n j g seen e q Ua | 60 — 15 — 15 = 30. 


4. The method 

We consider a (formal) differential equation 

V{z- F{z), F'(z), F"{z ),..., F^\z)) = 0, (4.1) 

where V is a polynomial with integer coefficients, which has a power series solution F{z) 
with integer coefficients. We assume that this power series solution F(z) is uniquely 
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determined by (14.1|) over the integers, and also over all rings (Z/p 7 Z), where p is a fixed 
prime number and 7 any positive integer. 

I 11 the above situation, we propose the following algorithmic approach to determining 
the series F(z) modulo a p-power pP a , for some non-negative integer a. We make the 
Ansatz 

p a + !~1 

m= E ai(z)$ l (z) modulo (4.2) 

i= 0 

with $( 2 ;) as given in (II.ip . and where the a^z)’s are (at this point) undetermined 
Laurent polynomials in z. Now we substitute (j4. 2[) into (14.ip . and we shall gradually 
determine approximations a it p(z) to a*(^) such that (14.11) holds modulo pP, for f3 = 
1,2,... ,p Q . To start the procedure, we consider the differential equation (14.11) modulo 
p, with 

P a+1 - 1 

F(z) = Fi(z) = ^ ai t i(z)^(z) modulo p. (4.3) 

i=0 

Using the elementary fact that <&'(z) = 1 modulo p, we see that the left-hand side of 
(14.ip is a polynomial in $(;?) with coefficients that are Laurent polynomials in z. We 
reduce powers $ fc (^) with k > p a+1 using the relation (which is implied by the minimal 
polynomial for the modulus p given in Proposition [£])□ 

(<L P (^) — <L(z) + z) P =0 modulo p p “. (4.4) 

Since, at this point, we are only interested in finding a solution to (14.ip modulo p, the 
above relation simplifies to 

<L pa+1 (^) — $ p “(^) + z pCl =0 modulo p. (4.5) 

Now we compare coefficients of powers <f> fc (z), k — 0,1,... ,p a+1 — 1 (see Remark HDD. 
This yields a system of p a+1 (differential) equations (modulo p) for the unknown Laurent 
polynomials a^i(z), % — 0 , 1 ,... ,p " +1 — 1 , which may or may not have a solution. 

Provided we have already found Laurent polynomials a^^z), i = 0,1,... ,p “ +1 — 1, 
for some (3 with 1 < /3 < p" — 1 , such that 

p a+1 - 1 

Y a iA z Wi z ) ( 2 * 4 * - 6 ) 

i=0 

2 Actually, if we would like to obtain an optimal result, we should use the relation implied by a 

minimal polynomial for the modulus in the sense of Section [2] But since we have no general 
formula available for such a minimal polynomial (cf. Item (2) of Remark [5] in that section), and since 
we wish to prove results for arbitrary moduli, choosing instead powers of a minimal polynomial for 
the modulus p is the best compromise. In principle, it may happen that there exists a polynomial in 

^(z) with coefficients that are Laurent polynomials in z, which is identical with F(z) after reduction 
of its coefficients modulo p p , but the Ansatz (14.211 combined with the reduction (14.41) fails because it is 
too restrictive. We are not aware of a concrete example where this obstruction occurs. The subgroup 
numbers of SL 2 CZ) which are treated modulo 8 in [T2J Section 11 ] by the method described here 

(specialised to p = 2), and modulo 16 by an enhancement of the method outlined in [12] Appendix D] 

are a potential candidate when considered modulo 2^ for /3 > 5. On the other hand, once we are 
successful using this (potentially problematic) Ansatz, then the result can easily be converted into an 
optimal one by further reducing the polynomial thus obtained, using the relation implied by a minimal 
polynomial for the modulus p p . 
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solves (14.1|) modulo p 13 , we put 

Oi,p+i(z) := aiA z ) + P 0 b it p +1 (z), i = 0,1, • • • ,P a+1 ~ 1, (4.7) 

where the fej i( g +1 (^)’s are (at this point) undetermined Laurent polynomials in Next 
we substitute 

P a+1 - i 

J2 a itP+1 (z)^(z) (4.8) 

i =0 

for F(z) in (14.ip . Using the fact that <3>'(z) = J2n=oP™ zP "* _1 m °dulo pP +1 , we expand 
the left-hand side as a polynomial in <h(^) (with coefficients being Laurent polynomials 
in z), we apply again the reduction using relation (14.41) . we compare coefficients of 
powers $ fc (^), k = 0,l,...,p" +1 — 1 (again, see Remark [TUP, and, as a result, we 
obtain a system of p a+1 (differential) equations (modulo jU +1 ) for the unknown Laurent 
polynomials bi t p+i(z), % — 0,1,... ,p a+1 — 1, which may or may not have a solution. If we 
manage to push this procedure through until (3 = p a — 1, then, setting a,i(z) = ai iP o(z), 
i — 0,1,... ,p a+1 — 1, the right-hand side of (14. 2|) is a solution to (14.11) modulo jf a , as 
required. 

Remark 10. As the reader will have noticed, each comparison of coefficients of powers of 
<f>(^) is based on the “hope” that, if a polynomial in <&(z) is zero modulo a p-power p 13 (as 
a formal Laurent series), then already all coefficients of powers of $(^) in this polynomial 
vanish modulo pr. However, this implication is false in general (see Lemma 39 in [ 12(1 
Appendix D] for the case of modulus 16). It may thus happen that the method described 
in this section fails to find a solution modulo p@ to a given differential equation in the 
form of a polynomial in <f>(z) with coefficients that are Laurent polynomials in z over 
the integers, while such a solution does in fact exist. As a matter of fact, this situation 
already occurred earlier in jl2 , Theorem 28] in the context of the analysis modulo 16 of 
the subgroup numbers of SX 2 (Z). There exists an enhancement of the method, outlined 
for the case of p = 2 in [T2], Appendix D], which (at least in principle) allows us to decide 
whether or not a solution modulo a given power in terms of a polynomial in <&(z) with 
coefficients that are Laurent polynomials in z over the integers exists, and, if so, to 
explicitly find such a solution. Since we do not make use of this enhancement in the 
present article, we refrain from presenting it here. 

It is not difficult to see that performing the iterative step (14.7ft amounts to solving 
a system of linear differential equations in the unknown functions bi t p + \(z) modulo p, 
where all of them are Laurent polynomials in z, and where only derivatives up to order 
p— 1 of the 6g3 + i(z)’s occur. Solving such a system is equivalent to solving an ordinary 
system of linear equations. This equivalence was explained for the case of p = 2 in 
m Lemma 12]. Since we do not need the general version of this lemma for arbitrary 
primes p in the present article, we leave it to the reader to work out the straightforward 
details. 

We remark that the idea of the method that we have described in this section has 
certainly further potential. For example, the fact that the series <&(z) remains invariant 
under the substitution z —» z p (or, more generally, under the substitution z —> z p , 
where h is some positive integer) — up to a simple additive correction — can be 
exploited in order to extend the range of applicability of our method to equations 
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where we not only allow differentiation but also this kind of substitution. This idea is 
actually already used in Section [71 see Theorem [HO 


5. The number of non-crossing connected graphs modulo 3-powers 


In this section, we apply our method to determining congruences modulo powers of 3 
for the numbers of non-crossing connected graphs with a given number of vertices. The 
origin of this analysis lies in a conjecture of Deutsch and Sagan [5] Conj. 5.18]. See 
(15.9j) and the subsequent paragraph for more information. 

A non-crossing graph on n vertices is a graph whose vertices are labelled by 1 , 2 , ,n, 
which can be embedded in the plane without crossings of edges in such a manner that 
the vertices are placed around a circle in clock-wise order. Let N n denote the number 
of non-crossing connected graphs with n vertices. It is known from [9] Eq. (16)] that 


N„ 


2 2n ~ 1 /§n - 2\ 
n \n — 1 J 



(5.1) 


Let N(z) := H n >i N n z n be the corresponding generating function. According to [6j 
Eq. (47)] (see also [8, Sec. 2.1]), the generating function N(z) satisfies the polynomial 
equation 

N 3 (z) + N 2 (z ) - 3z N(z) + 2z 2 = 0. (5.2) 

It is not difficult to see that (15.21) determines N(z) uniquely as a formal power series 
over Z, and as well modulo any power of 3, provided N 0 — 0 and N l — 1. 


Theorem 11. Let $(2) = ^ n>0 and let a be a non-negative integer. Then the 
generating function N(z), when reduced modulo 3 3<1 , can be expressed as a polynomial 
in <h(^) of degree at most 3“ +1 — 1, with coefficients that are Laurent polynomials in z 
over the integers. 


Remark. Computer experiments seem to indicate that the generating function N(z) in 
Theorem [TTJ when reduced modulo 3 3 °, can in fact be expressed as a polynomial in 
<L(z) of degree at most 3“ +1 — 1, with coefficients that are polynomials in z over the 
integers. 

Proof of Theorem [HI We apply the method from Section [H We start by substituting 
the Ansatz (j4.3[) in (15.21) and reducing the result modulo 3. In this way, we obtain 


3 “ +1 —1 3 a+1 -l 

<i0)$ 3i 0) + ^ <i0)$ 2 *0) 

2=0 2=0 

— a j,i(^) a i,i(^) < ^* +:, (^) + 2z 2 = 0 modulo 3. (5.3) 

0<j<?'<3 Q+1 -l 

We claim that the following choices solve the above congruence: 

00,1(2) = Sa(z) + s a (z ) modulo 3, 

a z a ,i( z ) = 1 — s a (z) modulo 3, 

02.3“, 1(2) = 1 modulo 3, 


(5.4) 
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where s a (z) = o^) with all other di } i(z) vanishing. In order to verify the claim, 

we substitute our choices in (I5.3H . For the left-hand side, we obtain 

$ 2 ' 3 “ +1 (4 + (i - 4(4) $ 3 “ +1 (4 + 4(4 + 4(4 

+ <h 4 ' 3 “(^) + (1 + 4(4 + s a (z )) $ 2 ' 3 “(4 + 4 ( 2 :) - 4 ( 2 :) + s 2 a (z) 

- (i - s a (z)) d> 3 ' 3 “(4 - (4(4 + «a(4) $ 2 ' 3Q (4 - M4 - 4(4) $ 3 “(4 + 2 4 

modulo 3. 

Using the relation (14.5p with p — 3 and reducing the obtained expression modulo 3, we 
arrive at 

<f> 2 - 3 “(4 - 24> 3 “(4 + 4' 3 “ + ( Sa (4 - 4(4) ($ 3 » - z 3a ) 

+ $ 3 » (<h 3Q (z) - 4“) + $ 2 - 3 » + 4(4 + 4(4 + 4(4 
- M4-4(4)$ 3 “(4 + 24 

= z 2 ' 3 " - z 3 a s a (z) + 2: 3 “4(4 + 4(4 + 4(4 + 4(4 - 4 modulo 3. 

By using the relation 

4(4 = s a(4 + 4 — modulo 3 (5.5) 

several times, this expression can be turned into 

z 2 ' 3 " - z 3a s a (z) + z 3 " ( s a (z ) + z 3a - 4 + (s«(4 + 4“ - 4 2 

+ s a (z ) ( s a (z ) + 4“ — 4 + 4(4 — 4 modulo 3. 

After expansion and reduction modulo 3, one sees that this expression reduces to zero. 

After we have completed the “base step,” we now proceed with the iterative steps 
described in Section [fl We consider the Ansatz (I4.6IMI4.8I) . where the coefficients aj i( g(4 

are supposed to provide a solution Fp(z) = j4=o 1 a i,p( z )^( z ) to (EZD modulo 3 /3 . 
This Ansatz, substituted in (15.2(1 . produces the congruence 

3“ +1 -i 3“+i_i 

- 3/3 E E %,444^+i(4* <+i (4 + F jK z ) + 444 - 3^44 + 22: 2 = 0 

j =0 i=0 

modulo 3^ +1 . (5.6) 

Since the sum has the prefactor 3^, we may reduce the 0^(4 modulo 3. By construc¬ 
tion, we have 

a 0 ,44 = a o,i(4 = 4(a) + s a (4 modulo 3, 

« 3 a ,p(z) = 0 , 3 °,i(z) = 1 — s a (z ) modulo 3, 

a 2-3“,/3(4 = a 2-3“,l(4 = 1 modulo 3, 
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and = 0 modulo 3 for all other z’s. If we substitute this in (I5.6j) and subsequently 
use (14.5|1 with p = 3 to reduce high powers of $(z), we obtain 

3“ +1 —1 3 Q+1 -1 

-3^(4(z) + s a (2)) Y Kp+i{ z W(?) s a {z)) Y kv+ 1 (-)$ 3 “ + *(» 

Z—0 Z—0 

3“ +1 -l 

-3* Y Kp+i(z)^ 2 ' 3a+i (z) + Fl(z) + Fl(z)-3zFp(z) + 2z 2 

i= 0 

3“ —1 

= ~ 3 ' 3 (( s «(a> + s «(a>) 

i =0 

- (1 - s Q (z)) z 3 °b i+2 . 3 “,/ 3 +i(z) - z 3a b i+3a ^ +1 (z)) $ z (z) 

2'3“—1 

- 3/3 Y (i S a( Z ) + Kp+l( z ) + (! - «a(2)) bi- 3ajP+1 (z) 

i= 3“ 

+ (1 — S a (,z)) bi +3 afi + i(z) + bi t / 3 + i(z) — z 3 &i+3“,/3+l(^)) ^ 1 { Z ) 

3 a+1 ~l 

-3^ Y {{s 2 a [z) + s a (z))b i)P+1 (z) + {I - s a {z))bi_ 3a ^ +l (z) 

i= 2-3“ 

+&i-2-3“,/3+l(^) + bi } /3+i(z)) <l>*(z) 

+ Fp(z) + Fp(z) — 3 zFp(z) + 2 z 2 = 0 modulo 3^ +1 . 

By our assumption on Fp(z), we may divide by 3 /3 . Comparison of powers of $( 2 ) then 
yields a system of congruences of the form 

M ■ b = c modulo 3, (5.7) 

where b is the column vector of unknowns (b t g + \ ( z ))i= o,i„. 5 3 a+i_! , c is a (known) column 

vector of Laurent polynomials in z, and M is the matrix 

(D(s 2 a (z) + s a (z)) D(-z 3a ) D(-z 3a ( 1 - 8 a (z))) \ 

B(l-S a (z)) D(s 2 a (z) + s a (z) + 1) D(1 - s a (z) - z 3a ) , 

\ D( 1) D(l-s a (z)) D(s 2 a (z) + s a (z) +1)/ 

with D(x) denoting the 3“ x 3" diagonal matrix whose diagonal entries equal x. We 
claim that 

det(M) = z 2 ' 3a modulo 3. 

In order to see this, we subtract s^(z) + s a (z ) times row i + 2 • 3“ from row i, i = 
0,1,..., 3“ — 1, and we subtract 1 — s Q (z) times row i + 2 • 3" from row i + 3 a , i — 
0,1,..., 3" — 1. If, in addition, we move rows 2 ■ 3“, 2 • 3“ + 1,..., 3 Q+1 — 1 to the top 
of the matrix, then it attains an upper triangular form, from which one infers that 

det(M) = (-1) 2 ' 3Q ' 3Q det 2 (D(s 3 q (z) - s a (z) - z 3 “)) det (D( 1)) 

= (4(*)-s«(*)-* 3 T 3a 

= z 2 ' 3a modulo 3, 


by means of (15.51) . as we claimed. As a consequence, the system (15.7j) is (uniquely) 
solvable. Thus, we have proved that, for an arbitrary non-negative integer a, the 
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algorithm of Section H] will produce a solution F 3 a(z ) to (15.21) modulo 3 3 ° which is a 
polynomial in $(2) with coefficients that are Laurent polynomials in z. □ 

We have implemented this algorithm. As an illustration, the next theorem contains 
the result for the modulus 27. 

Theorem 12. Let <I>(z) = J2 n > o^ 3 ”- Then we have 

Y N n z n = 18^ 3 + z 2 + z + (18z 3 + 12 z 2 )$(z) + (3z 2 + 15z)® 2 (z) 

n> 1 

+ (9 z 2 + 5z + 13)<L 3 (A) + (9z 2 + 6z + 24)<h 4 (^) + (15 z + 6)$ 5 (z) 

+ (18z + 4)$ 6 (^) + (18z + 21)$ 7 (z) + 12$ 8 (^) modulo 27. (5.8) 
hi [3 Conj. 5.18], Deutsch and Sagan conjectured that 

{ 1 (mod 3) if n = 3* or n = 2 • 3* for an integer i > 0, 

2 (mod 3) if n — 3* 1 + 3* 2 for integers A > A > 0, (5.9) 

0 (mod 3) otherwise. 

This conjecture was proved by Eu, Liu and Yell [7] by carefully analysing a binomial 
sum formula for N n modulo 3. A simpler proof was given by Gessel in [9j. It should be 
noted that his proof is essentially the one resulting from the a = 0 case of Theorem [Til 
Moreover, our results allow us to generalise this — in principle — to any power of 3. 
For the sake of illustration, we display the results for the moduli 9 = 3 2 and 27 = 3 3 . 

Corollary 13. The numbers N n of connected non-crossing graphs obey the following 
congruences modulo 27: 


(i) 


= 1 (mod 27) if. 

and only if n = 

1 , 

or n = 2-3* with 

i > 0 ; 


(ii) 

iVn 

= 2 (mod 27) i/ ; 

and only if n = 

10 

■ 3* with i > 1; 



hi) 

aa 

= 3 (mod 27) if. 

and only if, n = 

: 2 

• 3* + 1 with i > 1 

; or n 

= 3* + 2 with 


i > 

1, or n = 2 • 3 * 1 

+ 2 • 3 * 2 with A 

> 

A > 0 , or n = 4 

• 3 * 1 + 

3 * 2 + 3 * 3 with 


A - 

- 2 > A — 1 > A > 0 , or n = 3 * 1 + 

4- 

3*2 _|_ 3*3 A — 

3 > A 

-1 > A > 0, 


or n = 3 * 1 + 3 * 2 + 4 • 

3 * 3 with A — 3 > 

A 

- 2 > A > 0 ; 



IV) 

AA 

= 4 (mod 27) «/, 

and only if, n = 

3; 




(v) 

AA 

= 5 (mod 27) if, 

and only if, n = 

4- 

3* with i > 1; 



Vi) 

AA 

= 6 (mod 27) if, 

and only if, n = 

= 22-3* with i > 0, 

or n = 

3 H + 3*2 + ! 


with i\ — 1 > A > 2 , or n = 4 • 3 n + 3* 2 with A — 1 > A > 1, or n = 3 n + 4 • 3* 2 
with A — 2 > i 2 > 1, or n = 2 • 3 n + 3* 2 + 3* 3 with A — 2 > i 2 — 1 > A > 0, or 
n = 3 n + 2 • 3* 2 + 3* 3 with A — 1 > i 2 — 1 > i 3 > 0 ; or n = 3 n + 3* 2 + 2 • 3* 3 with 
A - 1 > A > A > 0 ; 

(vii) N n = 9 (mod 27) if and only if, n = 2 • 3 n + 2 • 3 * 2 + 3 * 3 with A > A > A > 0 , 
or n = 2 • 3 n + 3 * 2 + 2 • 3 * 3 with A > A > *3 > 0 , or n = 3 n + 2 • 3 * 2 + 2 • 3 * 3 
with A > A > A > 0 , or n = 2 • 3 n + 2 • 3 * 2 + 2 • 3 * 3 with A > A > A > 0 , or 
n = 3 n + 3 * 2 + 3 * 3 + 3 * 4 + 3 * 5 with A > A > A > A > A > 0 , 

or n = 2 • 3 n + 3 * 2 + 3 * 3 + 3 * 4 + 3 * 5 with A > A > A > A > A > 0 , 

orn = 3 n + 2 ■ 3 * 2 + 3 * 3 + 3 * 4 + 3 * 5 with A > A > A > A > A > 0 , or 

n = 3 n + 3 * 2 + 2 • 3 * 3 + 3 * 4 + 3 * 5 with A > A > A > A > A > 0 , or n = 3 n + 3 * 2 + 

3 * 3 + 2 ■ 3 * 4 + 3 * 5 withii > A > A > A > A > 0 , orn = 3* 1 +3* 2 + 3* 3 + 3* 4 + 2 -3 * 5 

with A > A > A > A > A > 0 ; 
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(viii) N n = 11 (mod 27) if and only if,n = 3* + 1 with i > 3; 

(ix) N n = 12 (mod 27) if, and only if, n — 7; or n = 40 • 3* with i > 0, or 

n = 2 • 3* 1 + 3* 2 + 1 with i\ — 1 > i 2 > 1, or n = 3 n + 2 • 3* 2 + 1 with 

11 — 1 > z 2 > 1, or n = 3 n + 3* 2 + 3* 3 + 3* 4 with i\ — 3 > i 2 — 2 > i 3 — 1 > i 4 > 0; 

(x) 7V n = 13 (mod 27) if and on// */, n = 3* with i > 2; 

(xi) 7V n = 15 (mod 27) f/, and only if, n — 3* + 4 with i > 3, or n = 4 ■ 3* + 1 with 

i >2, or n — 13 ■ 3* with i > 1, or n = 4 • 3 n + 2 • 3* 2 with i\ — 1 > *2 > 0, or 

n = 2 • 3 n + 4 • 3* 2 with R — 2 > i 2 > 0; 

(xii) 7V n = 18 (mod 27) if, and only if n = 2 ■ 3* 1 + 3* 2 + 3* 3 + 3* 4 with i 1 > 

1 2 > *3 > U > 0, or n = 3* 1 + 2 • 3* 2 + 3* 3 + 3* 4 with ii > i 2 > iz > U > 0, or 

n = 3 n + 3* 2 + 2 • 3* 3 + 3* 4 with ii > i 2 > 13 > ii > 0, or n = 3 n + 3* 2 + 3* 3 + 2 ■ 3* 4 
with ii > i 2 > *3 > *4 > 0, or n = 2 • 3 n + 2 • 3* 2 + 3* 3 + 3* 4 with i\ > i 2 > 

h > U > 0, or 7 i = 2 ■ 3 W + 3* 2 + 2 • 3* 3 + 3* 4 with ii > i 2 > i?, > 14 > 0, or 

n = 2-3 n +3* 2 +3* 3 +2-3* 4 withii > i 2 > h > U > 0, orn = 3 n +2-3 i2 +2-3 ?3 +3* 4 
with i\ > i 2 > ?3 > 74 > 0, or n = 3 n + 2 • 3* 2 + 3* 3 + 2 • 3* 4 with i\ > i 2 > 

1 3 > i\ > 0, or 7 i = 3 W + 3 * 2 + 2 • 3 * 3 + 2 • 3 * 4 with ii > i 2 > is > 14 > 0, or 

n = 3 n + 3 * 2 + 3 * 3 + 3 * 4 + 3 * 5 + 3 * 6 with i 1 > i 2 > i 3 > i 4 > i 5 > i 6 > 0; 

(xiii) N n = 20 (mod 27) if, and only if, n = 10, or 71 = 3 n + 3* 2 with R — 2 > i 2 > 1; 

(xiv) N n = 21 (mod 27) if, and only if, n = 7-3* with i > 1, or n = 4 ■ 3* 1 + 4 • 3* 2 

with ii — 2 > i 2 > 0, or n = 13 • 3* 1 + 3* 2 with ii — 1 > i 2 > 0, or 71 = 3 n +13 • 3* 2 
with i\ — 3 > > 0; 

(xv) N n = 23 (mod 27) if, and only if, n = 4; 

(xvi) N n = 24 (mod 27) if, and only if, n = 13, or n — 16 • 3* with i > 0, or 
n = 7 ■ 3 11 + 3 12 with i\ — 1 > i 2 > 0; 

(xvii) in the cases not covered by items (i)-(xvi), N n is divisible by 27; in particular, 
N n ± 7, 8,10,14,16,17,19, 22, 25, 26 (mod 27) for all n. 

Proof. We convert the right-hand side of (j5.8j) into a linear combination of series 
H ai ,a 2 ,...,a r {z) with all af s relatively prime to p , by means of Equation (13. 1 1) and Propo¬ 
sition El The result is 


N n Z n = 9 Z 2 + 15 z + (I82 + 13 )ddi(^) + ( 18 ^ + l)dd 2 (^) + ( 9 ^ + 20 )idij( 2 ) 

n> 1 

+ 127/2,1(2) + 127/12(2) + 247/1,1,1(2) + 127/4(2) + 37/2,2(2) 


+ 6/72,1,1(2) + 6/71,2,1(2) + 67/1,1,2(2) + 127/1,1,1,1(2) + 9/74,1(2) + 9/71,4(2) 

+ 9/72,2,1(2) + 9/72,1,2(2) + 97/1,2,2(2) + 187/2,1,1,1(2) + 187/1,2,1,1(2) 

+ 187/1,1,2,1(2) + 187/1,1,1,2(2) + 97/1,1,1,1,1(2) + 97/4,2(2) + 97/2,4(2) 

+ 187/4,1,1(2) + 187/1,4,1(2) + 187/1,1,4(2) + 97/2,2,2(2) + 187/2,2,1,1(2) + 187/2,1,2,1(2) 
+ 187/2,1,1,2(2) + 187/1,2,2,1(2) + 187/1,2,1,2(2) + 187/1,1,2,2(2) 

+ 97/2,1,1,1,1(2) + 97/1,2,1,1,1(2) + 97/1,1,2,1,1(2) 


Coefficient extraction following the algorithm described in Remark [H] then yields the 
claimed congruences. □ 


Considered modulo 9, the above corollary reduces to the following. 
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Corollary 14. The numbers N n of connected non-crossing graphs obey the following 
congruences modulo 9: 

(i) N n = 1 (mod 9) if, and only if, n — 1, or n — 2 ■ 3* with i > 0; 

(ii) N n = 2 (mod 9) if, and only if, n = 3* 1 + 3 12 with ii — 1 > i 2 > 0; 

(iii) N n = 3 (mod 9) if, and only if, n = 2-3 n +3* 2 withii > i 2 > 0, orn = 3 n +2-3* 2 
with i\ > i 2 > 0, orn = 2• 3 n + 2• 3* 2 with i\ > i 2 > 0, orn = 3 n + 3* 2 + 3* 3 + 3* 4 
with i\ > i 2 > *3 > ii > 0; 

(iv) TV™ = 4 (mod 9) */, and only if, n = 3* with i > 0; 

(v) iV n = 5 (mod 9) if, and only if, n — 4 • 3* with i > 0; 

(vi) N n = 6 (mod 9) if, and only if, n = 3 n + 3* 2 + 3* 3 with A > i 2 > i?, > 0, or 

n = 2 • 3 n + 3* 2 + 3* 3 with i\ > > ^3 > 0, or n = 3 n + 2 • 3* 2 + 3* 3 with 

ii > i 2 > > 0, or n = 3 n + 3* 2 + 2 • 3* 3 with i 1 > i 2 > i?, > 0; 

(vi) in the cases not covered by items (i)-(vii), N n is divisible by 9: in particular, 
N n ± 7, 8 (mod 9) for all n. 


Remark. In the first proof of the conjecture (15.9p in [7], Eu, Liu and Yeh worked from 
the sum expression 


N n = 


n 


2n—A 

E 

i=n —2 


3n — 3 
n + i + 1 J \n — 2 


which is provided in [8j Theorem 2], In order to achieve the proof, they considered four 
auxiliary sums of very similar nature, called there fi(n), i = 1,2, 3,4, of which they 
determined their behaviour modulo 3. These sums, and the numbers N n of connected 
non-crossing graphs, are reconsidered by Gessel in [9], who provides simpler proofs of 
the congruences modulo 3 established in [7], and also closed form expressions for the 
sums fi(n), i = 1,2,3,4. Gessel adds another sum, called / 5 (n), of which he also 
determines its behaviour modulo 3. 

All these sums can be treated in the same way as the numbers N n in Theorems [lTJ 
[T2l and Corollaries [131 El as we are going to explain now. In particular, in this way the 
behaviour of any of ff n), i — 1,2, 3,4, 5, modulo any power of 3 can be determined. 

The general family of sums introduced by Gessel in [9] is 


2n-\-j—k 

hj,k,i(n) = 

i=n—l 


f 3 n+j \ 
\n + i + k J 



? 


of which the sums fi(n), i — 1, 2, 3,4, 5, are the cases (j, k, l) = (1,1, 0), (0,1, 0), (0, 0, 0), 
(—1,1,1), (0,1,1), in this order. Gessel P, Sec. 3] proves that the generating function 

OO 

Hj,k,i(z) = 'y ^ hj t k,i(n)z 

n=k—j—l 


for these sums is given by 


(1 — 2 a(z)) 1 a k j l (z) 

(1 — 6 a(z) + 6 a 2 (z)) (1 — a(z )) k ~ 1 ’ 


where a(z) is the compositional inverse of ^(1 — z)( 1 -*2 z). The series a(z) thus being 
algebraic of degree 3, also all series Hj^jfz) are algebraic of degree 3. Below we list 
the algebraic equations satisfied by the generating functions for ffn), i = 1,2, 3,4, 5, 
and in each case we provide the base solution iq (^) for the modulus 3 3 “ (cf. (14.31) with 
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p = 3) for starting our method described in Section 01 The iterative steps carry through 
as well, as in the proof of Theorem [TIJ We leave the details to the reader. 

The generating function for the numbers fi(n) satisfies the equation 

(1 - 108z?)Hl lfi (z) - 3 H hlfi (z) + 2 = 0. (5.10) 

A base solution for the modulus 3 3a is Fi(z) = 1. The generating function for the 
numbers / 2 (n) satisfies the equation 

(1 — 108z 2 )hfQ 10 (z) — (1 + 9z)Ho } i t o(z) + z = 0, (5-11) 

and a base solution is Fi(z) = <h 3a (^) + s a (z). Next, the generating function for the 
numbers / 3 (n) satishes the equation 

(1 - 108 z 2 )Hl 0t0 (z) - (1 + 9z)H 0fii0 {z) -z = 0. (5.12) 

Here, a base solution is Fi(z) = — $ 3 “(^) — s a (z) + 1. The generating function for the 
numbers / 4 (rt) satishes the equation 

(1 - 108 z 2 )H* w (z) + (1 - 108 z 2 )H\ ltl (z) + 3z(l - 12z)H_ lxl (z) - 4z 2 = 0. (5.13) 

Since, modulo 3, this is the same equation as (15. 2[) . this generating function has the 
same base solution as N(z), namely the one given by (15.411 . Finally, the generating 
function for the numbers / 5 (n) satishes the equation 

(1 - 108 z 2 )Hl 1A (z) - Hq i i(z) - 8 z = 0. (5.14) 

Since, modulo 3, this is the same equation as (15.lip , (essentially) the same base solution 
(the difference occurring in the constant term), namely Fi(z) = <h 3a (^) + s a (z) + 1, 
applies. 


6. The number of Kreweras walks modulo 3-powers 

The subject of this section are congruences modulo powers of 3 for the numbers of 
so-called “Kreweras walks.” In his thesis p3j, Kreweras had considered (among many 
other things) a three-candidate ballot problem, where each of the three candidates, A, 
B, C, say, receives n votes. The problem which is posed is how many ways there are to 
count the votes so that at each point in time candidate A has at least as many votes as 
the number of votes for each of the other candidates. It is not difficult to see that this 
problem can be translated into the following problem of counting walks in the quarter 
plane: how many lattice walks in the plane from the origin to itself are there which 
consist of 3n steps from the set {(1,1), (—1,0), (0, —1)} and stay in the non-negative 
quadrant? Let K n denote the number of these walks. By definition, we set K 0 = 1. 
Kreweras [T3] proved that 

K„ = -^- -f 3 "j; (6.1) 

(n + l)(2n + 1) V n ) 

see also [T6j. To prove this formula turned out to be highly non-trivial, as well as the 
enumerative analysis of the more general question of how many walks there are which 
end at an arbitrary point in the non-negative quadrant. Kreweras [L4j and Niederhausen 
[T5 j obtained partial results in this direction. Finally, Bousquet-Melou [2] succeeded to 
compute the complete generating function for these numbers of walks. 
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Let K(z) := 'Yh n - > QK n z n be the generating function for the “Kreweras numbers” in 
(16.11) . From [2} Theorem 1], it can be derived that 

64s 2 K-\z) + 16s K 2 (z) - (72* - 1) K(z) + 54s - 1 = 0. (6.2) 

Clearly, this equation determines K(z) uniquely as a formal power series over Z, and 
as well modulo any power of 3. 

Theorem 15. Let <F(s) = an d let a be a non-negative integer. Then the 

generating function K(z), when reduced modulo 3 3a , can be expressed as a polynomial 
in ^(s 1 / 2 ) of degree at most 3" +1 — 1, with coefficients that are Laurent polynomials in 
s 1//2 over the integers. 

Proof. We apply the method from Section [4] We start by substituting the Ansatz (14. 3 j) 
with z replaced by s 1//2 in (j6.2[) and reducing the result modulo 3. In this way, we 
obtain 


/Y. 1 al(z^mz^) + z 3 Y 

i=0 i=0 

-z Y, «L 1 (^ 1/2 )%, 1 (^ 1/2 )^ +i (^ 1/2 ) 

0<«<j<3“+ 1 -l 

3 a+1 -l 

+ — 1 = 0 modulo 3. (6.3) 

i=0 

We claim that the following choices solve the above congruence: 

a 0 ,i(s 1/2 ) = s _1 s 2 (s 1 / 2 ) modulo 3, 

a3 a ,i(^ 1//2 ) = — s _1 s Q ,(s 1 ^ 2 ) modulo 3, 
a2-3“,i(^ 1 ^ 2 ) = z ~ l modulo 3, 

where s a (s) is the same polynomial as in Section 0 with all other aj ;1 (s 1//2 ) vanishing. 
In order to verify the claim, we substitute our choices in (16.31) . For the left-hand side, 
we obtain 

s 2 (s- 3 <L 2 ' 3a+ V /2 ) - s- 3 s 3 (s 1/2 )$ 3 “ +1 (s 1/2 ) + W 3 4(s 1/2 )) 

+ s (s” 2< h 4 ' 3a (s 1//2 ) + z-^Kz^^iz 1 / 2 ) + z-^ 2 )) 

- s (-s- 2 ^^ 1 / 2 )^ 1 ^ 1 / 2 ) + s" 2 s 2 (s 1 / 2 )^^ 1 / 2 ) - z- 2 s 3 a (z l / 2 )&‘*(z 1 / 2 )') 

+ z-^^iz 1 / 2 ) - z-'s^z 1 / 2 )^^' 2 ) + z-'sKz 1 ' 2 ) - 1 modulo 3. 
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Using the relation (14. 5 p with p — 3 and z replaced by z 172 , and reducing the obtained 
expression modulo 3, we arrive at 

(z-^v 72 ) + + 

- (z~hl(z 1/2 )$ 3a (z 1/2 ) - z 3a/2 -h 3 a (z 1/2 )) +Z~ 1 S 6 a ( Z 1/2 ) 

+ (z-^'^iz 1 / 2 ) ~ Z/ 2 -^ 3 ^^/ 2 )) + 2 ) + Z^S^ 1 / 2 ) 

+ (z^z 172 )* 3 ^ 2 ) -z 3 */ 2 - 1 ^^/ 2 )) 

- z^s 2 ^ 172 )* 2 ^^ 172 ) + z-^l^ 2 )^^ 2 ) 

+ z _1 <h 2 ' 3 “(z 1/2 ) - z _ 1 s q (z 1 / 2 )$ 3 “(z 1/2 ) + z _1 s 2 (z 1/2 ) - 1 modulo 3. 

By collecting terms, this expression simplifies to 

z 3 “—1 + z 3 “ /2_1 s 3 (z 1/2 ) + z-^z 172 ) + z-^z 172 ) 

- z 3 “ 72 - 1 s a (z 172 ) + z-^ 2 (z 172 ) - 1 modulo 3. 

By repeatedly using the relation (15.51) with z replaced by z 172 , this expression can be 
turned into 

z 3 “-i + z 3 " 72 - 1 (s Q (z 172 ) + z 3 “ 72 - z 172 ) + 3-1 ( Sa {z 1 ' 2 ) + z 3 “ 72 - z 172 ) 2 

+Z- 1 Sa (z 172 ) ( Sq (z 172 ) + Z 3 “ 72 - ^/2)_.3V2-1 Sq( M/2 ) + .- 1 s 2 ( M/ 2 ) _ 1 modulo 3 . 

After expansion and reduction modulo 3, one sees that this expression reduces to zero. 

After we have completed the “base step,” we now proceed with the iterative steps 
described in Section [H We consider the Ansatz (j4.6p ~ fj4.8p with z replaced by z 172 , 
where the coefficients a^p{z 172 ) are supposed to provide a solution 

3«+! _1 

F p{z) = Y a^{z 1/2 )^\z 112 ) 

i =0 

to (15.21) modulo 3^. This Ansatz, substituted in (15.2j) . produces the congruence 


3<2+i_ 1 3 Q +i_ 1 3 Q: +i_l 

-3»z y , y v w*- i/2 )*v i/2 ) 

j =0 i=0 i =0 

+ 64z 2 F|(z) + 16zF|(z) — (72z — l)Fp(z) + 54z — 1 = 0 modulo 3^ +1 . (6.4) 

Since the sum has the prefactor 3^, we may reduce the aj i( a(z 172 ) modulo 3. By con¬ 
struction, we have 

ao i/3 (z 172 ) = a 0i i(z 1//2 ) = z _1 s 2 (z 1/2 ) modulo 3, 
a 3 “,/?(z 172 ) = a 3 “,i(z 172 ) = —z _1 s Q (z 172 ) modulo 3, 

<3'2-3“, / s(^ 1//2 ) = 02 - 3 “,i(z 172 ) = z' 1 modulo 3, 
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and = 0 modulo 3 for all other z’s. If we substitute this in (16.411 and subsequently 
use (I4.5p with p = 3 and z replaced by z 1 / 2 to reduce high powers of ^(z 1 / 2 ), we obtain 


2—0 2—0 

3«+i_X 3«+i_2 

-3 » J2 &< J , +1 L 1/2 )'i>' +2 - 3 "(A 2 ) + 3S V & iij8+1 (W 2 )4><(A 2 ) 


i=0 i=0 

~2 773/^3 , 172 / 


3 “ —1 


i=0 


+ 64z 2 F|(z) + 16zFp(z) — (72z — 1 )Fp{z) + 54z — 1 
-3” Y, (sl(z 1/2 )b t ,?Mz' /2 ) + 2 3 " /2 2>4* I/2 )'' i+ 2.3.,W,(* 1/2 ) 

-2 s “ /2 t i+ 3.,j + i(.- 1/2 ) - w* ,/2 )) ms 12 ) 

-3® Y (s 2 ( 2|1/2 ) (, i,/>+l( 2l/2 ) “ »o(^ 1/2 )ii-3“J>+l(2 1/2 ) 

-S„(3 1/2 )6 j+3 ., /i+1 (3 1/2 ) - 3 3 " /2 6i + 3»,J + l(^ 1/2 )) 4'(i 1/2 ) 


2 ' 3 “ —1 


i=3 a 


3 Q+1 -1 


- 3/3 X] (s 2 a (z 1/2 )b i>f)+1 (z 1/2 ) - s a (z 1/2 )bi_ 3a!f}+1 (z 1/2 ) + 6i_ 2 .3“,/3+i(- 1/2 )) $\z 1/2 ) 


i=2-3° 


+ 64z 2 F|(z) + 16 zFp(z) — (72z — l)Fp(z) + 54 z — 1 = 0 modulo 3 /i+1 . 

By our assumption on Fp(z), we may divide by 3A Comparison of powers of <L(z 1//2 ) 
then yields a system of congruences of the form 


M ■ b = c modulo 3, (6.5) 

where b is the column vector of unknowns (6,: :/ 9 +i ( z .^a+i-!, c is a (known) 

column vector of Laurent polynomials in z, and M is the matrix 

(D {si (z 1 / 2 ) - 1) D(-z 3a / 2 ) D(z 3 a / 2 s a (z 1 / 2 )) \ 

D (- S „(A 2 )) D ( s 2 a (z V 2 )) D i-sUz 1 ' 2 ) - z 2 "/ 2 ) , 

V 0(1) o(- s „(W 2 )) O( s 2 (0/ 2 )) j 

with D(x) denoting the 3 a x 3“ diagonal matrix whose diagonal entries equal x, as 
before. In the same manner as in the proof of Theorem [lTJ on sees that 

det(M) = z 3a modulo 3. 

As a consequence, the system (16. 5p is (uniquely) solvable. Thus, we have proved that, for 
an arbitrary non-negative integer a, the algorithm of Section U] will produce a solution 
F 3 a(z) to (15.2p modulo 3 3 ° which is a polynomial in ^(z 1 / 2 ) with coefficients that are 
Laurent polynomials in z. □ 


We have implemented this algorithm. As an illustration, the next theorem contains 
the result for the modulus 27. 


Theorem 16. Let <f>(z) = J2 n>0 z 3n . Then we have 
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K n z n = 1 + 12 <S> 2 (z 1/2 ) + 11 z- l,2 $\z 1/2 ) + 6^ _1 $ 4 (^ 1/2 ) + 15^" 1/2 $ 5 (^ 1/2 ) 

n> 0 

+ ^" 1 $ 6 (^ 1/2 ) + ?>z- l <$>\z 1/2 ) modulo 27. (6.6) 

Corollary 17. The numbers K n of Kreweras walks obey the following congruences 
modulo 27: 

(i) K n = 1 (mod 27) «/, and only if, n = 0; 

(ii) K n = 2 (mod 27) if, and only if, n — 1; 

(iii) Ji n = 3 (mod 27) if, and only if, n = 3* with i > 1, or n = |(13 • 3 ? — 1) with 

i > 2, or n = |(3* +11) with i > 4, orn = 2 • 3* +1 with i > 3, orn = 20 • 3* — 1 

withi > 1, or n = |(3 n +3* 2 + 3* 3 + 3* 4 — 2) withi\ — 3 > i 2 — 2 > i 3 — 1 > i 4 > 1; 

(iv) Ji n = 5 (mod 27) «/, and only if, n = |(3 n + 3* 2 — 2) with i\ — 2 > i 2 > 2; 

(v) Ji n = 6 (mod 27) if, and only if, n — 3* 4- 1 with i > 3, or n = |(3 n + 3* 2 ) 

with i\ — 2 > i 2 > 1, or n — |(7 • 3* 1 + 3* 2 — 2) with ii — 1 > i 2 > 1; or 

n = |(3 n + 7 • 3* 2 — 2) with i\ — 1 > i 2 > 1; 

(vi) Jv„ = 7 (mod 27) «/, and on/y if, n — 3* — 1 with i >2; 

(vii) Jv„ = 8 (mod 27) if, and only if, n — A, or n = 2 • 3* — 1 with i > 2; 

(viii) Ji n = 9 (mod 27) if and only if, n = 3 n + 3* 2 + 3* 3 — 1 with i\ > i 2 > *3 > 1, 

or n = |(2 • 3 W + 3* 2 + 3* 3 + 3* 4 + 3* 5 — 2) with i\ > i 2 > i 3 > ii > i$ > 1, or 

n = |(3 n + 2 • 3* 2 + 3* 3 + 3 M + 3* 5 — 2) with i\ > i 2 > i 3 > *4 > i*> > 1, or 

n = |(3 n + 3* 2 + 2 • 3* 3 + 3* 4 + 3* 5 — 2) with *i > *2 > *3 > *4 > *5 > 1, or 

n = |(3 n + 3 12 + 3* 3 + 2 • 3* 4 + 3* 5 — 2) with i\ > i 2 > i 3 > u > *5 > 1, or 

n = |(3 n + 3* 2 + 3* 3 + 3* 4 + 2 • 3* 5 — 2) with i\ > i 2 > is > ii > i$ > 1; 

(ix) K n = 12 (mod 27) 2/, and only if n — ^(4 • 3* 1 + 3* 2 — 1) with i\ — 1 > i 2 > 2, 

or n = |(3 n + 4 • 3* 2 — 1) irit/z ii — 2 > i 2 > 2, or n = |(3 W + 3* 2 + 2) 

with ii — 2 > i 2 > 3, or n = |(13 • 3* 1 + 3* 2 — 2) with i\ — 1 > i 2 > 1, or 

n — 2 • 3 n + 2 • 3* 2 — 1 with ii — 2 > i 2 > l, or n = |(3 n + 13 • 3* 2 — 2) with 

ii - 3 > i 2 > 1; 

(x) Jv„ = 14 (mod 27) if, and only if, n = |(3* + 1) with i > 4, or n = 5 • 3* — 1 
with i > 2; 

(xi) K n = 15 (mod 27) if, and only if, n — |(7 • 3* — 1) with i >2, or n = |(3* + 5) 

with i > 2, or n — 2 ■ 3* wn£/i i > 1, or n = 11 • 3* — 1 with i > l, or 

n = i(2-3 n +3* 2 +3* 3 —2) withii—2 > i 2 —1 > i 3 > 1, orn = |(3 n +2-3* 2 +3* 3 —2) 
with ii — 1 > i 2 — 1 > i 3 > 1, or n = |(3 n + 3* 2 + 2 • 3* 3 — 2) with i\ — 1 > i 2 > 

* 3 > i; 

(xii) A'„ = 16 (mod 27) if, and only if, n = 2; 

(xiii) = 17 (mod 27) if, and only if, n = 5; 

(xiv) K n = 18 (mod 27) if, and only if or n = |(2 • 3* 1 + 2 • 3* 2 + 3* 3 + 3* 4 — 2) 
with ii > i 2 > i 3 > i 4 > 1, or n = |(2 • 3* 1 + 3* 2 + 2 • 3* 3 + 3* 4 — 2) with 

h> i 2 > h> U> 1, or n = |(2 • 3 n + 3* 2 + 3* 3 + 2 • 3* 4 — 2) with i x > 

12 > is > U > 1, or n = |(3 n + 2 • 3* 2 + 2 • 3* 3 + 3* 4 — 2) with ii > i 2 > i 3 > 

ii > 1, or n = i(3 n + 2 • 3* 2 + 3* 3 + 2 • 3* 4 — 2) with i\ > i 2 > i 3 > i 4 > 1, 

or n = |(3 n + 3* 2 + 2 • 3* 3 + 2 • 3* 4 — 2) with ii > i 2 > i 3 > i 4 > 1, or 

n = |(3 n + 3* 2 + 3* 3 + 3* 4 + 3* 5 + 3* 6 — 2) with i\ > i 2 > i 3 > i 4 > i 5 > i 6 > 1; 
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(xv) K n = 21 (mod 27) if and only if, n = 19, or n = 3 n +3* 2 —1 with ii > i 2 > 1, or 

n = i(3 ?1 +3* 2 +3* 3 —1) with i\ —2 > i 2 — 1 > i 3 > 2, orn = -f(4-3 n +3* 2 +3* 3 —2) 

with ii — 2 > i 2 — 1 > is > 1, or n = |(3 n + 4 • 3* 2 + 3* 3 — 2) with — 3 > 
*2 — 1 > *3 > 1, or n = -(3* 1 + 3* 2 + 4 • 3* 3 — 2) with R — 3 > i 2 — 2 > i 3 > 1; 

(xvi) K n = 23 (mod 27) if, and only if, n = 14; 

(xvii) K n = 24 (mod 27) «/, and only if, n = 10, or n = |(2 • 3* 1 + 3* 2 — 1) with 
i 1 — l>i 2 >2, orn = |(3 n + 2 -3* 2 — 1) with i\> i 2 > 2, or n = 3 n + 2 - 3* 2 — 1 
with i\ — 2 > i 2 > 1, or n = 2 • 3 n + 3* 2 — 1 with R > i 2 > 1; 

(xviii) K n = 26 (mod 27) «/, and only if, n = ^(3* — 1) with i > 3; 

(xix) m the cases not covered by items (i)-(xvii), K n is divisible by 27; in particular, 
K n ^ 4,10,11,13,19, 20, 22, 25 (mod 27) for all n. 

Proof. By means of Equation (13. ip and Proposition [9l we convert the right-hand side 
of (16.6p into a linear combination of series 77ai,a 2 ,...,ar (^ 2 ) with all afs relatively prime 
to p. The result is 

Y,K n z n = 21z~ 1/2 H 1 {z 1/2 ) + 7z~ x H 2 {z xl2 ) + 5 z~'H ltl {z x/2 ) 

n> 0 

+ 9z~ 1/2 H 2tl (z 1/2 ) + 9z~ 1/2 H lt2 (z 1/2 ) + 18z- 1/2 H 1Aj1 (z 1/2 ) 

+ 3 z- l H A (z l/2 ) + 21 Z - 1 H 2 , 2 {z 1/2 ) + 15 z- l H 2 ^(z 1/2 ) 

+ lbz-'Hwiz 1 ' 2 ) + 150- 1 77i i i, 2 (^ 1/2 ) + 3^- 1 i^i,i,i,i(^ 1/2 ) 

+ 9z~ 1 H i2 {z 1 ^ 2 ) + 9z- 1 H 2A (z l/2 ) + 18^- 1 i^ 4 ,i,i(^ 1/2 ) + 18 z' 1 H lA1 (z l/2 ) 

+ 18z“ 1 i7i i i )4 (2; 1/2 ) + 9z _1 i7 2)2>2 (^ 1 ' /2 ) + 18z~ x H 2t2AA (z 1/2 ) 

+ 18 z- l H 2X2A {z l/2 ) + 18 z-'H^^z 1 ' 2 ) + 18z _1 i7 1)2)2)1 (,2^ 2 ) 

+ 18£ -1 i/ li2i i i2 (z 1/2 ) + I8z _1 H lAfia {z 1/2 ) + 9^ _1 iy 2 ,i,i, 1 ,i(^ 1/2 ) 

+ 9^ _1 i^i, 2 ,i,i,i(^ 1/2 ) + 9z~ 1 H lyl)2AA (z 1/2 ) + 9z~ 1 H 1Xh2A (z 1/2 ) 

+ 90 _1 77i i i i i.i i2 (^ 1 ' /2 ) + lSz -1 #!^^!^!^ 1 / 2 ) modulo 27. 

Coefficient extraction following the algorithm described in Remark [S] then yields the 
claimed congruences. □ 

If we restrict Corollary [T7] to modulus 3, then it reduces to the following simple 
assertion. 

Corollary 18. The number K n of Kreweras walks of length 3 n is congruent to 1 mod¬ 
ulo 3 if, and only if, n — 3* — 1 with i > 1, it is congruent to 2 modulo 3 if, and only 
if, n — |(3 n + 3* 2 — 2) with i\ > i 2 > 0, and it is divisible by 3 in all other cases. 

7 . FUSS-CATALAN NUMBERS MODULO p-POWERS 

In this section, given a positive integer h and a prime number p , we determine the 
behaviour of the Fufi-Catalan numbers 

F(»i *):=!(*") (7.1) 

n \n — 1 J 
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modulo powers of p, provided k is itself a power of p, say, k = These numbers have 
numerous combinatorial interpretations; cf. [TJ pp. 59-60]. 

By using the Lagrange inversion formula (see [2U, Theorem 5.4.2]), it is easy to 

see that the generating function f p -h(z) = 1 + X^n>i n(n-i) zn satisfies the functional 
equation 

z f^ h {z) - f P -,h( z ) + 1 = 0. (7.2) 

It is straightforward to verify that this equation has a unique formal power series solu¬ 
tion over Z, and as well over any power of p. 

In order to determine the coefficients of f p; h( z ) modulo powers of p, we have to use 
a variant of the series < L(^) defined in (II.ip . namely 

= X> P " V(ph-1) - (7-3) 

n> 0 

The theorem below generalises Theorem 33 in [12] to arbitrary prime numbers. 

Theorem 19. For a prime number p and a positive integer h, let $ p -h( z ) be the series 
defined in (17.31) , and let a be a further positive integer. Then the generating function 
fp-h(z), when reduced modulo pP a , can be expressed as a polynomial in <& p ^{z) of degree 
at most p( a+1 ) /l — with coefficients that are Laurent polynomials in z l ^ p _1 - > over the 
integers. 

Proof. For ease of notation, we replace z by z ph in (17.21) . thereby obtaining the equa¬ 
tion 

~ UA Z ) + 1 = 0, (7.4) 

with f p -h(z ) = f p -h(z ph ~ 1 ). We now have to prove that, modulo p p °‘ h , the series f v -h{z) 
can be expressed as a polynomial in 

OO 

*RfcW = £ ~-" h (7.5) 

n= 0 

of degree at most pl a+1 ') h — l 5 with coefficients that are Laurent polynomials in 0. 

It is readily verified that 



^ P p . h (z) - %, h {z) + z = 0 modulo p, 

(7.6) 

whence 

( f I+l(^) - %-,h(z) + .?) =0 modulo f ah . 

(7.7) 

We modify our 

Ansatz (|4.2[) to 



p(a+l)h_Y 



fpA z )= a X z WpA z ) modulo p p °‘ h , 

(7.8) 


i=0 


where the afizf s are (at this point) undetermined Laurent polynomials in z. 

3 In principle, one could use the generalisations of Lucas’ theorem due to Davis and Webb [4 ], 
and to Granville m, respectively, to analyse the explicit expression for the Fufi-Catalan numbers 
modulo a given p-power. However, this approach would be rather cumbersome in comparison with our 
method, and it is doubtful that one would be able to derive a result on the same level of generality as 
Theorems [19] [701 or Corollaries [7TJ [77] 
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Next, we gradually find approximations a^p{z) to a,i(z) such that (I7.4j) holds modulo 
pP, for /3 = 1,2,... ,p ah . To start the procedure, we consider the functional equation 
(17.4p modulo p, with 

p(a + l)h_^ 

fp-,h(z)= a iA z )K -, h ( z ) modulo p. (7.9) 

i= 0 

It is readily verihed that the choice of 

ao,i (*) = 

k =0 

a pa h A (z) = z _1 , (7-10) 

with all other aj.i(z)’s equal to zero indeed leads to a solution of (j7.4[i modulo p. 

After we have completed the “base step,” we now proceed with the iterative steps 
described in Section [4] Our Ansatz here (replacing the corresponding one in (14.611 — (I4.8[) ) 
is 

p(a + l)h _ 1 

fp-,h(z) = a hP+i( z )®l-,h( z ) modulo pP +1 , (7.11) 

1=0 

with 

ai,p+i(z) := di,p(z) + p\p +l (z), i = 0,1,... ,p (a+1)ft - 1, (7.12) 

where the coefficients a,i,p{z) are supposed to provide a solution 

hi*) = a iA z )%,h( z ) 

2=0 

to (17.41) modulo p&. This Ansatz, substituted in (17.41) . produces the congruence 

p(a+l)h_i 

z ph ~ l ff{z) -fp(z) +pp WWW) + 1= 0 modulo pP +1 . (7.13) 

2 = 0 

By our assumption on fp(z), we may divide by pP. Comparison of powers of $ P] h(z) 
then yields a system of congruences of the form 

bi,p+i(z) + Poli(^) = 0 modulo p, i — 0,1,... — 1 ; (7.14) 

where Poh(^), i = 0,1,... i p^ a+1 ) h — l 5 are certain Laurent polynomials with integer 
coefficients. This system being trivially uniquely solvable, we have proved that, for an 
arbitrary positive integer a, the modified algorithm that we have presented here will 
produce a solution f pa h(z) to (17.41) modulo p p ° which is a polynomial in <3> p -h(z ) with 
coefficients that are Laurent polynomials in z. □ 

It should be observed that the a = 0 case of the above proof (see in particular 
(17.101) 1 shows that F(n\p h ) = 1 (mod p) for n = ( p hl — l)/(p — 1), i = 0,1,..., and 
F(n, p h ) = 0 (mod p) otherwise. 

We have implemented the algorithm contained in the above proof. As an illustration, 
we display below the result obtained for h — 1 and the modulus p 2 . This result was 
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first guessed from the automatically obtained results for p = 3,5, 7, but, once found, it 
is easily verified directly by substitution in (17.21b 

Theorem 20. Let 3>(z) = Y2 n >o zP " ■ Then we have 


l) "" = P$ P-1 (2 1/(P_1) ) ~(P~ l)^" 1/(p " 1) $ p (^ 1/(p - 1) ) 

n> 1 ' ' 7 

+ pz - 1/ (p- 1) $2 P - i ( z i/( P - i) ) modulo p 2 . (7.15) 
(Human) inspection reveals that this expression can actually be drastically simplified. 
Corollary 21. Let <E>(z) = ■ Then we have 

~ (j™ = z -V(p- 1 )$p( z 1 /(^- !)) modulo p 2 . (7.16) 

n>l ' ' 7 


Proof. One applies relation (j4.4j) with a = 0 and z replaced by ^ to the terms in 
(17.15j) which have coefficient p. □ 

Explicitly, this leads to the following congruences modulo p 2 . 

Corollary 22. The Fufi-Catalan numbers F(n;p ) = 7 J obey the following congru¬ 
ences modulo p 2 : 

(i) If n = —C (p l — 1) with i > 1, then F(n;p) = 1 (mod p 2 ). 

(ii) V 

[aip 11 + a 2 p 12 4-b a r p lr - l) 


n = 


p — 1 

with all ai ’s relatively prime to p, a\ + a 2 + ■ ■ ■ + a r = p, r > 2, and ii > i 2 > 
• • • > i r > 0, £/ien we /iaue 


F(n; p) = 


p! 


ai! 02 ! • ■ ■ a r ! 


(mod p 2 ). 


(iii) In i/ie cases not covered by items (i) and (ii), the Fufi-Catalan number F(n\p) 
is divisible by p 2 . 

Proof. By means of Equation (13.ip and Proposition [91 we convert the right-hand side 
of (15.8p into a linear combination of series i7 ai ,a 2 ,...,a r , with all af s relatively 

prime to p. The result is 


r( pn 1) * n = - _1/(p ' 1) ^i(- 1/(p ' 1) ) -1 

n>l'7 / 


+ 2 


-l/(p-l) 


E E 

r=l ai,...,a r >l 

aiH-|-a r =p and r>2 


p! 


ai! 02 ! • • • a r ! 


77ai,a 2 ,...,a r (^ 1/(P 1} ) modulo p 2 . 


Coefficient extraction then leads to our claim. 


□ 
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8. The number of blossom trees modulo p-powers 


The combinatorial objects which we treat in this section are so-called “blossom trees.” 
These are a particular kind of trees which are of great significance in the combinatorial 
understanding of the enumeration of maps; see e.g. the survey [3]. The particular 
blossom trees that we are interested in are the ones in [[18j Sec. 3]. Since the precise 
definition is slightly technical and not needed here, we omit it, and instead refer the 
reader to [fl8j . 

For an odd positive integer k, let B{n\ k ) be the number of blossom trees constructed 
from k-axy trees with n white nodes by adding a black node with k —1 buds on each inner 
edge. Let B k (z) := _|_ ^2 n>l B(n\ k) z n be the corresponding generating function. 

Schaeffer [T8j Cor. 2] proved that 


B(n ; k) 


k + 1 f kn \ 
n((k — l)n + 2) \n — lj ' 


( 8 . 1 ) 


As we show in the Appendix, the generating function B k (z) 
equation 


satisfies the polynomial 


(fc+l)/2 


£ (-1) S 77- 


k + 1 


s =0 


k — s + 1 


(k — s + l)(k — s) 

i)‘ f— i— V 'sfcW + f-T 


,k—2s+l | 

(k-1 


{ 2 

k + 1 i 

(k- T 

2 ' 

< 2 




k -1 


= 0 . ( 8 . 2 ) 


It is not difficult to see that (I8.2jl determines B k (z) uniquely as a formal power series 
over Z, and as well over any power of a prime number p as long as k ^ 1 (mod p). 


From now on, let p be a fixed odd prime number. 


Theorem 23. Let $(x;) = an d let a be a non-negative integer. Then the 

generating function B p (z), when reduced modulo pP a , can be expressed as a polynomial 
in $ (z 1 /( p-1 ' ) ) of degree at most p a+1 — 1, with coefficients that are Laurent polynomials 
in over the integers. 


Proof. Again, we apply the method from Section [4l To begin with, we need a “base 
solution” 

P a+1 - i 

F 1 (z)= (8-3) 

i =0 

to (18.21) modulo p. At this point, it is useful to observe that, using Fermat’s little 
theorem, the functional equation (j8.2j) . when taken modulo p, reduces to 

z 2 B p p (z ) - 2 ^ +1 l/ 2 zB^ +1 ^ 2 {z) + B p {z) - 2" 1 = 0 modulo p, (8.4) 

where here, and in the following, 2 _1 denotes the inverse of 2 modulo p. We claim that 
the following choices in (18.3ft solve the above congruence: 

a 0A {z) = 2 - 1 z~ 2/(p - 1) sl{z 1/(p - 1) ) modulo p, 


a pa l (z) = z 2 ^ p 1 - ) s^z 1 ^ 1 )) modulo p, 
d 2 -p a ,i(z) = ir\- 2 ^ modulo p , 
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where s a (z) = z pk , with all other a,i t i(z) vanishing. In order to verify this claim, 

we hrst observe that, with the above choices, F\{z) can be expressed as a square, namely 

F,{z) = 2 -^- 2 /h-D ($p“(^V(p-i)) + Sa (^/(P-D)) 2 . 

We substitute this in (18.4[) . For the left-hand side, we obtain 

2 -U-2/( P -i) ($p“(^i/(p-D) +Sq (.i/(p-D)) 2p 

_ Z - V ( P-D ($p“( 2 i/(p-D) + Sq (^/(p-D))P +1 

+ 2-W- 2 /(p-!) ($P Q (^/(P-1)) + ^(^/(p- 1 ))) 2 - 2 _1 modulo p. 

Using the relation (14. 5 p with z replaced by and reducing the obtained expression 

modulo p, we arrive at 

2 -i z -2/(p-i) ^$p“ +1 ( z i/(p-i)) +sP(^/(p-i)^ 2 

_ z -2/(p-l) ^$P Q+1 ^1/(P-1)) ($P“( Z 1/(P-1)) +5 a (z 1 /( p - 1 >)) 

+ 2- 1 W 2 /(p- 1 ) ($p q (^/(p-D) + Sq (^/(p-D)) 2 _ 2-1 
= 2-U- 2 /(P-!) - ^P a /(P-1) + s Q ,_|_ 1 (^!/(P—!)) - ^ 1 /(P- 1 )) 2 

- z- 2 ^-V (r^- 1 )) - + Sa+1 (^ 1 /( p - 1 )) - ^/Cp-D) 

x ($p“( 2 1/( p- 1) ) +Sq ( 2 1/( p- 1) )) 

+ 2- 1 ^- 2 /(p- 1 ) ($p q (^/(p-D) + Sa ( z i/(P-D)) 2 - 2- 1 
= 2- 1 z" 2/(p - 1) ($P“(z 1/(p - 1} ) + s a (z l/{ - p - l) ) - z 1/{p - l) ) 2 

- Z ~ 2 ^~!) - ^ 1/(p - 1} ) 

x ($P“(^ 1/( P- 1) ) +Sa (^ 1/( P- 1) )) 

+ 2- 1 W 2 /( p - 1 ) (^“(^/(p- 1 )) + s q (z 1 Z(p- 1 ))) 2 - 2' 1 modulo p. 

By collecting terms, this expression may be simplified to zero. 

After we have completed the “base step,” we now proceed with the iterative steps 
described in Section [4] Our Ansatz here (replacing the corresponding one in (I4.6H — (j4.8jl ) 
is 

p(o + l)/i „ 1 

B p (z) = ^ a it p + i(z)$ l modulo p^ +1 , (8.5) 

i =o 

with 

a%,p+i{z) := ^{z) +pPbij+\{z), i = 0,1,... ,p (a+1)h - 1 , (8.6) 

where the coefficients a^(z) are supposed to provide a solution 

p(o + l)h _ 1 

Ff, M= W (P^- 11 ) 

i=0 
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to (18.21) modulo p&. This Ansatz, substituted in (18.21) , produces the congruence 


-,\(p+l)/2 


2(-l) 


p — 1 


( P ~ l)/2 


P + 1 


n“ +1 -l 


d “ +1 -1 


(P-1)/2 


[ £ 6 i , s+1 (x)<I»U I/( '- I) ) ] [ V 

i=o 


p — 1 


i=0 

p-1 /p Q+1 -l 




t=0 


(P+l)/2 


++ E (-d't— 


P + 1 


s=0 


(p — s + l)(p — s) \ s 


p-s+1 y- 2s+ i 


p — i 




p —i 


Fb(- 


P + 1 f P ~ 1 


= 0 modulo pP +l . (8.7) 


Since the sum has the prefactor p&, we may reduce the a^ 0 {z) modulo p. By construc¬ 
tion, we have 

a O}0 (z) = a 0 ,i(z) = 2~ 1 z~ 2 / <J> ~ 1 ^Sa(z 1 ^- 1 ' ) ) modulo p, 
a P <*,p{z) = a p * t i(z) = z~ 2/<(P ~ l) s a {z 1/{ - p ~ l) ) modulo p, 

a 2 -p a ,p{ z ) — 02 -p“,i(^) = 2 - 1 Z ~ 2 ^ P ~ 1 ' 1 modulo p, 

and di^i = 0 modulo p for all other V s. Equivalently, we have 

F 0 {z) = 2~ 1 z~ 2 ^ (^(^/(p- 1 )) + s^z 1 ^)) 2 modulo p. 

If we substitute this in (18.71) and divide both sides of the congruence by pP (which is 
possible by our assumption on Fp(z)), we obtain 

- \ ” e ‘ W ^ b 170 " 1 ’)) (^"( AM ) + s „ p i/ &>- 1 )))-- 1 


i=0 


'P Q+1 -1 

E kj>+iM*V /< ^ I) ) ) + G f (z) 

i=0 


p a+1 -l 


- £ 6 i , ( , +1 (z)*U I/< " _I) ) E(-l) 3 « 3!, “( 


i =0 


/ p—1 

E< 

d=o 


z 1/(p - 1] ) s p - j - 1 {z l/{p ~ l) ) 


p“ +1 -l 


+ b i,0+i( z W{ zl/iP lj ) +G 0 {z) = o modulo p, 


i =0 


where Gp(z) is some polynomial in <f> (z 1 ^ p 1 - > ) whose coefficients are Laurent polyno¬ 
mials in ^/(p-d with integer coefficients. Upon using (I4.5j) with z replaced by z l ^ p ~ 1 ' > 
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to reduce high powers of ^(z 1 ^ we convert this congruence into 

- E' E~ (*)<&-' 

3= 0 i=jp a 

P-1 (j+l)P“-l 

-EE 4>v i/(p - 1) )(-i 

j —0 i=p a 
p-1 jp a -1 

+ E E «‘(‘ 1/<i, - 1) ) i+Mp »,„ + 1 ( Z K^- 1 (P /(p - 1) ) 

J=0 2—0 

( p«+l_l \ 

y, 6 ,: i/ 3 + i(^)<f> , (^ 1/(p_1) ) I + = 0 modulo p, 

Comparison of powers of < f > (^) then yields a system of congruences of the form 

M ■ b = c modulo p, ( 8 . 8 ) 


where b is the column vector of unknowns (2))i=0,l,.- ^a+Lj, c is a (known) column 

vector of Laurent polynomials in z, and M is the matrix 


(D(l - s p a - 1 ) 

D{z a ) 

D(-z A s a ) 

D(z A s p a ~ 3 ) 

D( — Z A Sa 2 ) ^ 

D(s p ~ 2 ) 

Di-s^ 1 ) 

D(s a + z A ) 

... D{-s p - 3 - z A sP- A ) 

D{s p ( p 2 + z a s P cA 3 ) 

D(-s p a ~ 3 ) 

D(s p - 2 ) 

Di-sP- 1 ) 

... D(s p a ~ 4 + z A s p - 5 ) 

D(-s p a - 3 - z A s p a - 4 ) 

D(s a ) 

D(sl) 

D(s 3 a ) 

Di-s/sr 1 ) 

D(s a + z A ) 

V ^(-i) 

D(s a ) 

D(sl) 

D(s p a - 2 ) 

Di-sP- 1 ) / 


where s a is short for , A is short for p a /(p— 1), and D(x) denotes the p a xp a 

diagonal matrix whose diagonal entries equal x. More precisely, the (r, t)-block of the 
matrix, 0 < r, t < p — 1, is given by 

'Dil-sPp 1 ), if r = f = 0, 

D((—l) t ~ 1 z A s t ~ 1 ), if r = 0 and t > 0, 

D((— l) r_t+1 s^ _1 ~ r+i ), if r > t but not r — t — 0, 

v L)((-l) t_r_1 s^ _r_1 (s a + z A )), if 0 <r <t. 

We claim that 

det(M) = z pa modulo p. 

In order to see this, we add 1—s^ _1 (^) times row i+(p—l)-p a to row i,i — 0,1,... ,p a — 1, 
and we add times row i + p a to row i, i = p a ,p a + 1,... ,p a+1 — 1. This 

leads to the matrix 


( ° 

D(s a + z A - s p ) 

* 

* 

* \ 

0 

0 

D(s a + z A - s p ) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

D(s a + z A - s p a ) 

\D(- 1) 

D{s a ) 

D(-s 2 a ) 

••• D(s p ~ 2 ) 
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If, in addition, we move rows (p — 1) • p a , (p — 1) • p a + 1,..., p a+1 — 1 to the top of the 
matrix, then it attains an upper triangular form, from which one infers that 

det(M) = (-1 )(P-V-P a -P a det+ ^“/(p-i) _ s p (z 1/(p - 1} ))) det (£>(1)) 

= (saO 1 ^- 1 )) +^“ / (p- 1) - s p(z 1/ (p- 1 ))) (p_1) ‘ p “ 

= z pa modulo p, 

as was claimed. As a consequence, the system (j8.8[) is (uniquely) solvable. Thus, we 
have proved that, for an arbitrary non-negative integer a, the algorithm of Section [4] 
will produce a solution F pa (z) to (j8.2[) modulo p pa which is a polynomial in 
with coefficients that are Laurent polynomials in z. □ 

A direct corollary of the a = 0 case of the preceding proof is the following. 

Corollary 24. The blossom tree numbers B(n]p ) = n (( p ^t)n +2 ) (rf-i) °^ e ?/ the following 
congruences modulo p: 

(i) If n = ipzi(p n + p 12 — 2) with i\ — 1 > i 2 > 0, then B(n;p) = 1 (mod p). 

(ii) Ifn = ffi-j- (2 p l — 2) with i > 1, then B(n]p) = ^ (mod p). 

(iii) In the cases not covered by items (i) and (ii), the number B{n]p) is divisible 
by p. 

Obviously, one can do better by using an implementation of the algorithm contained 
in the proof of Theorem [23] As an illustration, we display below the result obtained for 
the modulus p 2 . Again, this result was first guessed from the automatically obtained 
results for p = 3,5, 7, but, once found, it is easily verified directly by substitution in 

4H2D- 

Theorem 25. Let <h(^) = J2 n >o zI>n ■ Then we have 


^l + J2 B (n;p) z' 

n> 1 


( P + l)z- 1/{p - 1 '>$(z 1/(p - 1 ')) - P±l z -2/(p- 


_l_ z - 2 /(p-i)^p+ir z i/(p-i)^ 


1 )$ 2 (^ 1 /( p_1) ) 

modulo p 2 . (8.9) 


Explicitly, this means the following. 


Corollary 26. The blossom tree numbers B(n\p) = n ((JfX) n + 2 ) (n^i) °^ e ?/ the following 
congruences modulo p 2 \ 

(i) If n = _ ((p -f i)p* — 2) with i > 0, then B{n\p ) = 1 (mod p 2 ). 

(ii) Ifn = ffi-j- (2 p l — 2) with i > 1, then B(n\p) = (mod p 2 ). 

(iii) If n = —A (p n + p 12 — 2) with i\ — 1 > i 2 > 0, then B(n;p ) = p + 1 (mod p 2 ). 

(iv) If 

n = ——- Uip n + a 2 p 12 H-b - 2) 

p — 1 

with all ai’s relatively prime to p, Gq + a 2 + • • • + a r — p + 1, r > 2, and 
ii > *2 > ■ ■ ■ > i r > 0, then we have 


B{n\p) 


(P+ 1)! 


ai! a 2 \ ■ ■ ■ a r \ 


(mod p 2 ). 
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(v) In the cases not covered by items (i)-(iv), the number B(n;p) is divisible by p 2 . 


Proof. By means of Equation (13 .1 p and Proposition [HI we convert the right-hand side 
of (18. 9p into a linear combination of series H ai , a 2 ,..., a r with all af s relatively 

prime to p. The result is 


^2 B (n-,p) 


z n = P±l z -2/(p-l) H 2 ( z l/(p-l)j + ( p+ i) z -V(p-i)H ljl (z 1 / ( r- 1] ) 


n > 0 


pz~ 2 ^ v ~ 1 " > H p+ i (^ 1 /( p_1 )) 


E 


P+1 

+ z -i/(.r-DJ2 

r =1 ai,...,a r >l 

cti H- \-a r =p -\-1 

r> 2 and a^p for all i 


- ^ — :H ai a 2 a (z 1 ^ p 1 - ) ) modulo 

a x ! a 2 ! • ■ ■ a r \ 11 2 ’"’’ A ’ 


p 2 . 


Coefficient extraction then yields the claimed congruences. 

Appendix: Proof of the functional equation (j8.2[) 
Let k be a positive integer and let 


□ 


nw = E - 

< ^ n 


n > 1 


kn 


n \n 


be the generating function for the general Fufi-Catalan numbers. We claim that the 
blossom tree generating function B k (z) can be expressed in terms of T k (z) as follows. 


Lemma 27. For all positive integers k, we have 


B k (z) = (1 +T k (z)) k 


k- 1 


(1 +T k (z)) 


(A.l) 


Proof. We verify that the coefficient of z n is the same on both sides of (1A.ll) . for 
n — 0,1,... Clearly, the constant coefficient equals (k + l)/2 on both sides. 

Now let n > 1. It is easily seen by Lagrange inversion (see [20, Theorem 5.4.2]), that 
the series T k (z) satisfies the equation 

T k (z) = z(l + T k (z)) k . (A.2) 

In other words, it is the compositional inverse of z/{ 1 + z) k . By applying Lagrange 
inversion again, we infer 

k- 1,„ . m / 1 , _ 1X , (1 + z) kn 


W) (1 + T t W) (k - L-T (1 + Tk (z)fj = I (z- 1 ) (1 - (fc - 1)*) 


z" 


= - ({/- 1 } (1 + z) kn -{k- 1) (E 1 ” 2 ) (1 + z ) kn ) 
1 ( ( kn \ ,, ( kn 

~{k~ 1) 


n 


kn 

n — 1 
k + 1 


n((k — 1 )n + 2) \n — 1 


kn 


n — 2, 

= B(n\ k). 


□ 
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Proof of (18.211 . In order to establish (18.211 . we substitute the right-hand side of (I A. 1 1) 
for Bk{z ) in (18.21) . After applying the binomial theorem, we obtain 


3=0 


1 I . | k k - j 


Jk-\ 


(fc+l)/2 s 

E E<-i) ,+) 

s=0 j =0 


X(z) k+j 
k + 1 


(k — s + l)(k — s) 
k — s + 1\ (s 


u _i \ s +t 




A; - 1 


fc-i 


fc-1 


*(*) ( * - —X(z) I + (-1) 2 


k + 1 / k — 1 


k -1 


on the left-hand side, where X(,s) = 1 + Tifz). In order to simplify, we use the relation 
(IA.2[) in the first term, we replace s by s—j, and then write the sum over j in standard 
hypergeometric notation 


F 

p-Tq 


Cl 1 5 • • • 5 &p m 


E 


(^l)r 


„y, m! (*l)r 


•p)m 

/C 


( 6 . 


q Jm 


where the Pochhammer symbol (a) m is dehned by (a) m := a(ct + 1) • • • (a + m — 1), 
m > 1, and (a) 0 := 1. In this manner, we arrive at 


E(-iy Q k k -> (X(z) - 1 )X(z)‘ 


+ 2(-l) 


k-\-l 


k- 1 


zX{z) k+1 + {-lfk(k + 1) 


ifc - 1 


fe-i 




fc-1 


1 y {k + l) (k - s - l)\ k _ s+1 fk-1 
^ 1 j (A; — 2s + 1)! s! 

■zX(z) s 3 F 2 


3—0 


k — s, —| + 1 
--s-L-^^-s + l’ 1 

9! '’WJ- 


L2 






The 3 F 2 -series can be evaluated by means of the Pfaff-Saalschiitz summation formula 
(see m (2.3.1.3), Appendix (III.2)]) 


a, b, —n 

c, 1 + a + b — c — rC 


(c - a) n (c - b) n 
(i c) n {c-a-b) n ' 


3 F 2 
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provided n is a non-negative integer. Thus, we obtain 

k ,, 


Eh) 

j=0 


j ( k ) k k-j 


fc- 1 


*(*(*) - l)X(z)' 


+ 2(-l) fc+1 ( zX(z) k+1 + (-l) k k(k + l) fk 1 


fc-i 


zX(z) 


fc -1 




( k + s - 1) /fc-1 


fc* 


~ (-1)* 
k 

EM ) 

j=o 


s=0 

fc-1 


s—1 


fc + 1 


zl(z) 6 


fc-1 


A(z) fc 


fc-1 


X(z) + (-1)' 


.fc + 1 /fc — 1 


fc -1 


./fc 


J M fc fc_j 


, / fc -1 


z{X{z) - 1)X(^ 


+ (-l) 


fc+i 


fc -1 


zX(z) k+1 + (-1)' 


,fc + 1 /fc — 1 


v—', . s (fc + S — 1). /fc — 1 

2J- 1 )'-—«— ~ k 


s=0 


2 

s—1 


2 

fc + 1 


fc-i 


zx( z y 


zX(z) 


-(-i y 


fc-i 


fc-i 


fc-i 


(fc-—A^)J +(-l) 2 


fc + 1 /fc — 1 


fc-i 


In the last expression, the sums can in fact be extended to run over all integers j 
respectively s. A comparison of powers of X(z) then makes it obvious that the sum 
over j cancels with the sum over s. To the remaining terms A" fc+1 (z) and X k (z) we 
apply again the relation flA.21) . This turns the above expression into 


- 1 ) 


k +1 


fc-i 


-iy 


{X{z) - l)X(z) + (-I)' 


.fc + 1 /fc — 1 


fc-i 


fc -1 


fc-i 


X(z) k 


fc -1 


X(z) +(-l) 


(X(z) - 1) 


.fc + 1 /fc — 1 


fc -1 


= 0 , 


which finishes the proof of ( 18.21) . 


□ 
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